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ABSTRACT 


THREE  RESONANT  NONLINEAR  OPTICAL  PROCESSES 
IN  ATOMIC  CESIUM 


by 

Arlee  Virgil  Smith 

Chairman:  John  F.  Ward 

We  have  investigated  three  resonant  nonlinear  pro- 
cesses in  atomic  cesium  vapor.  They  are  third  harmonic 
generation,  two  photon  absorption,  and  parametric  fre- 
quency conversion.  Resonant  enhancement  of  each  is 
provided  by  tuning  the  ruby  laser  near  half  the  6S  - 9D 
transition  frequency  of  28828.90  cm  1.  The  parametric 
process  involves,  in  addition,  a resonance  with  one  of 

the  P states. 

■ > 

The  experimental  apparatus  consists  of  a tuneable, 
narrow  bandwidth,  Q-switched,  ruby  laser  with  maximum 
power  of  1 MW,  a sample  cell  containing  .1  torr  or  less  of 
cesium  vapor,  and  beam  optics  and  optical  detectors. 

We  have  studied  the  variation  of  each  of  the  three 
processes  with  laser  intensity  and  frequency.  At  exact 
resonance  the  THG  is  found  to  be  proportional  to  the 
square  of  the  laser  intensity  and  the  TPA  is  proportional  - 


to  laser  intensity  over  a large  power  range.  The 
resonance  lineshape  for  both  processes  is  found  to  be 
power  dependent.  A detailed  analysis  of  these  two 
processes  based  on  steady  state  density  matrix  calcula- 
tions and  including  A.  C.  Stark  shifts,  photoionization, 
and  beam  geometry  is  shown  to  accurately  describe  most 
of  our  observations.  ^ 

Signals  have  been  observed  at  348  nm,  362  nm,  389  nm, 
459  nm,  3 ym  and  1.4  ym  corresponding  to  transitions  from 
10P,  9P,  8P,  and  7P  to  6S  and  from  9D^  to  7P  and  8P 
respectively  and  have  been  shown  to  be  produced  by  a four 
photon  parametric  process.  The  signals  are  radiated 
colinearly  with  the  laser  beam  and  the  four  high  frequency 
signals  attain  maximum  powers  of  a few  watts.  Exponential 
gain  has  been  measured  for  each  at  low  powers  and  com- 
pared with  calculated  values  of  gain.  Saturation  is 
observed  in  each  case.  We  show  that  this  is  probably 
caused  by  populating  the  P states.  The  effect  of  9D^ 
population  on  the  value  of  the  gain  coefficient  is  also 
discussed  and  shown  to  be  important. 


pip.  / (TflfF-T i-ei&t'i 


THREE  RESONANT  NONLINEAR  OPTICAL  PROCESSES 


IN  ATOMIC  CESIUM  * 


i r»  0 'Tfi  I s / - i 


Arlee  Virgil^ Smith 


A dissertation  submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy 
(Physics) 

in  The  University  of  Michigan 
1977 


Doctoral  Committee : 


Professor  John  F.  Ward,  Chairman 
Professor  C.  Wilbur  Peters 
Associate  Professor  Richard  L.  Sears 
Professor  Yukio  Tomozawa 
Associate  Professor  William  L.  Williams 


44*  '146 


I 

I 

I 

T 

i 

I 

I 

l 


| 

ACKNOWLEDGEMENTS 

I thank  Professor  John  Ward  for  providing  me  the 
opportunity  to  perform  this  research  and  for  suggesting 
the  topic.  His  physical  insight  and  patient  guidance 
throughout  have  been  invaluable. 

Kang  Leung's  help  in  the  early  stages  of  the  work 
is  gratefully  appreciated.  I thank  also  Irving  Bigio 
for  his  advice  and  example  while  an  officemate. 

A special  thanks  to  friend  and  colleague  Charles 
Miller  for  his  assistance  and  cooperation  throughout 
and  for  his  sense  of  humor. 

I am  especially  grateful  to  my  wife  Elaine  for 

cheerfully  contributing  many  hours  of  work  to  the  prep- 

* 

eration  of  this  manuscript  and  for  her  continued  support 
and  encouragement. 

Finally,  the  financial  support  furnished  by  the 
United  States  Air  Force  Office  of  Scientific  Research 
is  acknowledged.  wc.  ^ ^ 3-  -V  ftcoafc  _'_h~ 

1 uJttWWO 

1 JUS  ! '/■v,nN> 

I ^ 

k 

ii 


□ n' 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGEMENTS  ii 

LIST  OF  TABLES  V 

LIST  OF  FIGURES  vi 

LIST  OF  APPENDIXES  viii 

CHAPTER  1.  INTRODUCTION  1 

CHAPTER  2.  THE  EXPERIMENTAL  APPARATUS  8 

2.1  The  Cesium  Cell  8 

2.2  The  Attenuator  13 

2.3  The  Detector  & Monitor  Boxes  ...  13 

2.4  The  Data  Acquisition  System  ....  17 

2.5  The  Laser  17 

CHAPTER  3.  BASIC  THEORY  OF  TWO  PHOTON  ABSORPTION, 
THIRD  HARMONIC  GENERATION  AND  PARA- 
METRIC GAIN  26 

3.1  Two  Photon  Absorption  30 

3.2  Third  Harmonic  Generation  36 

3.3  Parametric  Gain  42 

3.4  Summary  63 

CHAPTER  4.  THIRD  HARMONIC  AND  TWO  PHOTON  ABSORP- 
TION EXPERIMENTS  65 

4.1  Steady  State  Saturation  Model  ..  67 


4.2  Comparison  of  Model  and 
Experimental  Results  . . 


CHAPTER  5. 


APPENDIXES 


4.3  Conclusions  

PARAMETRIC  FREQUENCY  CONVERSION  . . 

5.1  The  Data  

5.2  Comparison  of  Data  and  Theory 

5.3  Limitations  of  Theory  and 

Modifications  

5.4  Summary  and  Conclusions  


REFERENCES 


LIST  OF  TABLES 


Table  Page 

3.1  List  of  definitions  and  conventions  ....  28 

3.2  Matrix  elements  for  atomic  cesium  60 

3.3  Calculated  parametric  gain  coefficients 

and  threshold  powers  60 

5.1  Comparison  of  measured  and  calculated 

power  and  density  parametric  gain 
coefficients  and  threshold  powers  118 


v 


LIST  OF  FIGURES 


Figure  Page 

1.1  Energy  level  diagram  of  atomic  cesium  ..  6 

2.1  Experimental  layout  9 

2.2  Cesium  density  calibration  data  12 

2.3  Fluorescence  detection  system  16 

2.4  Laser  cavity  and  tuning  controls  19 

2.5  Frequency  dependence  of  laser  tuning 

elements  20 

2.6  Frequency  calibration  data 23 

3.1  Graph  of  TPA  cross  section  35 

3.2  Graph  of  THG  efficiency  41 

3.3  Graph  of  <X^>  for  xg  = 0 50 

3.4  Graph  of  <Xj>  for  xg  = 1 50 

3.5  Graph  of  <XS>  for  Xg  = 0 51 

3.6  Graph  of  <Xg>  f°r  x9  = 1 51 

3.7  Graph  of  <Xj><Xs>  for  xg  = 0 54 

3.8  Graph  of  <Xj><Xs>  for  xg  = 1 54 

3.9  Lineshape  of  parametric  gain  for  Xg  = 0.  56 

3.10  Lineshape  of  parametric  gain  for  xg  = 1.  56 

4.1  Measured  power  laws  for  TPA  and  THG  ....  66 

4.2  Calculated  lineshapes  for  TPA 72 

4.3  Calculated  lineshapes  for  THG 72 

4.4  Calculated  power  dependence  of  TPA  76 

4.5  Calculated  power  dependence  of  THG  77 


vi 


Figure  Page 

4.6  Power  law  of  THG  with  Stark  shift  79 

4.7  Oscilloscope  traces  of  fluorescence  ....  83 

4 . 8 Fluorescence  variation  with  cesium 

density  84 

4.9  Hyperfine  structure  of  atomic  cesium  ...  87 

4.10  TPA  and  THG  measured  lineshapes  (single 

mode  laser)  90 

4.11  Fit  of  computed  and  measured  TPA  line- 

shapes  (single  mode  laser)  92 

4.12  Fit  of  computed  and  measured  THG  line- 

shapes  (single  mode  laser)  95 

4.13  TPA  and  THG  measured  lineshapes  (multi- 

mode  laser)  97 

4.14  Measured  power  laws  of  TPA  and  THG  98 

4.15  Fit  of  computed  and  measured  line- 

shapes  (multimode  laser)  99 

4.16  Computed  TPA  power  law  101 

4.17  Computed  THG  power  law  102 

5.1  Observed  signals  and  idlers  112 


5.2  Temporal  behavior  of  parametric  signal  . 112 

5.3  Variation  of  10P  parametric  signal  with 


laser  power  114 

5.4  Variation  of  9P,  8P,  and  7P  parametric 

signals  with  laser  power  115 

5.5  Variation  of  8P  and  10P  parametric 

signals  with  cesium  density  116 

5.6  Lineshape  of  10P  parametric  signal  120 


r 


vii 


t 


CHAPTER  1 

INTRODUCTION 

The  electric  polarization  induced  in  matter  by  an 
applied  electric  field  is,  in  general,  nonlinear  in  the 
applied  field.  Consequently,  for  sufficiently  small 
electric  fields,  it  is  useful  to  express  the  polarization 
as  a power  series  in  the  field. 

p = y + xxE  + x2e2  + x/  + •••  (l.D 

Here  P represents  the  total  electric  polarization,  y 
the  permanent  polarization,  E the  applied  field,  and 
the  coefficients  (x's)  are  called  susceptibilities. 

The  nonlinearity  of  P has  many  interesting  conse- 
quences. For  example,  if  x2  is  nonzero  and  the  electric 

field  has  a cos (at)  time  dependence,  then  the  term  X2E 

, 2 
will  produce  polarization  proportional  to  cos  (at)  which 

has  Fourier  components  at  the  frequencies  2a  and  0 . 
Likewise,  the  fourth  term  in  the  expansion,  X3E^»  will 
produce  polarization  at  the  frequencies  3a  and  la,  and 
higher  order  terms  will  generate  various  other  harmonics. 
If  the  applied  field  has  more  than  a single  Fourier  com- 
ponent, for  example  E = E^cos(a1t)  + E2cos(a2t),  then 
the  term  X2E  produce  polarization  at  frequencies 

(a^  + a2)  and  (a^  — tu2 ) as  well  as  0,  2a^ , and  2a2-  The 


1 


2 


fourth  term  will  generate  polarization  at  (2oj^  + ) , 

(2co^  - ) » (2^2  + cu1 ) , and  so  on. 

The  polarization  which  is  induced  at  each  point  in 
the  nonlinear  material  in  turn  radiates  an  electromagnetic 
wave  with  the  same  set  of  Fourier  components  as  the  po- 
larization. The  study  of  this  phenomenon  when  the  applied 
field  has  optical  frequency  components  constitutes  the 
field  of  nonlinear  optics.  Besides  sum  and  difference 
frequency  generation,  such  effects  as  self -induced 
transparency,  optical  nutation,  Raman  scattering,  and 
multiphoton  absorption  are  also  included.  Reviews  of  the 
subject  may  be  found  in  articles  by  Giordmaine , Franken 
and  Ward,  and  Bloembergen,  among  others.1 

The  susceptibilities,  although  written  as  constants 
in  eqn.  1.1,  must,  in  general,  be  frequency  dependent 
tensors  in  order  to  adequately  relate  applied  electric 
fields  of  various  directions  of  polarization  and  fre- 
quencies to  the  induced  polarization  of  another  polari- 
zation and  frequency.  The  conventional  notation  for  X 
is  illustrated  in  eqn.  1.2. 

e“  - xijk<-»;»1.“2)  e7  Ek2  <l-2> 

The  subscripts  here  refer  to  the  direction  of  polariza- 
tion (i,j,k  = #,y,z)  and  the  superscripts  to  the  fre- 
quency. Each  combination  of  Fourier  components  of  the 


applied  field  has  a x associated  with  it. 


If  the  nonlinear  sample  material  is  an  atomic  gas, 
we  can  be  more  specific  about  the  properties  of  the  x's* 
The  symmetry  of  the  sample  requires  that  all  x ' s with  an 
odd  number  of  subscripts  be  zero  and  only  certain  tensor 
elements  are  allowed  to  be  nonzero  for  x ' s with  an  even 
number  of  subscripts.  To  cite  a familiar  example,  only 
diagonal  elements  of  the  linear  susceptibility  x^j  are 
nonzero  and  they  are  equal  (x  = x = x ) • 

j axx  Ayy  az  z 

The  form  of  the  frequency  response  is  also  known 
for  an  atomic  gas.  For  the  example  x^j  » the  response 

is  described  by  the  resonance  denominator  = (ft-cu) 

. i v . 

where  fi  is  oi  - Here  co  is  a natural  or  resonance 

o 2 o 

frequency  and  p is  a damping  factor.  For  X^j^/  the 
response  is  determined  by  the  denominator  D,  = (ft  - co.  ) 

j d I 

(ft,  - id-  - to.  ) (ft  - a),  - a)-  - w. ) .2  This  form  leads 
b 2 1 C j 2 1 

to  a variety  of  interestinq  phenomena  includ- 
ing resonantly  enhanced  sum  and  difference  frequency 
generation  processes,  two  photon  absorption  or  emission, 
and  stimulated  Raman  scattering. 

The  present  work  is  concerned  primarily  with  re- 
sonant processes  mediated  by  in  an  atomic  gas. 

Specifically,  we  have  observed  three  resonant  nonlinear 
processes  in  atomic  cesium  vapor.  At  least  two  of  the 
applied  fields  are  supplied  by  a ruby  laser  in  each  in- 
stance. The  atomic  cesium  - ruby  laser  system  is 


convenient  for  the  investigation  of  resonant  processes 
because  the  laser  is  tunable  over  a range  which  in- 


cludes the  frequency  (Egg  - EgD.  )/2fi,  where  Egg  and  EgD 
are  the  energies  of  the  6S  and  9D^  states  of  atomic 
cesium  (see  fig.  1.1).  This  means  that  the  second  term 
in  the  resonance  denominator,  D3 , can  be  made  very 
small.  Each  of  the  three  processes  and  its  relationship 
to  this  two  photon  resonance  is  described  briefly  be- 
low. 


The  first  process  is  two  photon  absorption. 


When  the  laser  is  tuned  to  half  the  6S  - 9D^  transi- 


tion frequency  of  cesium,  two  laser  photons  can  be  simul 
taneously  absorbed  by  a ground  state  (6S)  cesium  atom 
inducing  a transition  to  the  9D^  level.  This  nonlinear 


absorption  is  the  consequence  of  having  an  imaginary 


part  in  the  same  way  that  linear  absorption  is  caused 


by  the  imaginary  part  of  x • In  this  case  the  imaginary 


part  of  x^  is  large  at  resonance  because  the  middle 


term  in  the  resonant  denominator  Dg  is  imaginary. 

The  second  process  of  interest  is  third  harmonic 
generation.  In  this  case,  tuning  the  laser  near  the 
two  photon  resonance  provides  enhancement  in  the  pro- 
duction of  the  third  harmonic  of  the  laser.  Both  real 


and  imaginary  parts  of  x^  contribute  to  this  process. 


The  energy  levels  involved  in  third  harmonic  generation 
are  indicated  in  fig.  1.1. 
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The  third  process,  parametric  frequency  conversion, 
can  be  regarded  as  two  coupled  difference  frequency 
generating  processes.  One  mixes  the  laser  field  with 
itself  and  with  what  is  conventionally  called  the 
idler  field  to  produce  polarization  at  the  signal  fre- 
quency a)  where  cu  = 2to  - oj  (see  fig.  1.1).  The  other 

u b XI 

mixes  the  laser  field  with  itself  and  with  the  signal 
frequency  field  to  produce  the  idler  frequency  polariza- 
tion. Under  some  circumstances  the  coupling  between 
these  processes  can  provide  positive  feedback  so  that 
exponential  growth  of  the  signal  and  idler  waves  occurs. 
The  end  result  is  that  energy  is  exchanged  between  the 
laser  field  and  the  signal  and  idler  fields.  Two 
resonances  can  provide  enhancement  of  this  process.  One 
of  them  is  the  two  photon  resonance  as  before  but  the 
second  is  the  single  photon  6S  - nP  resonance,  where 
n,  the  principle  quantum  number,  can  be  10,  9,  8,  7,  or 
6 (see  fig.  1.1) . 

Each  of  these  three  nonlinear  processes  has  been 
previously  observed  in  a variety  of  nonlinear  materials 
including  atomic  vapors.*  In  particular,  Leung  has  ob- 
served each  in  cesium  vapor  using  much  the  same  experi- 
mental apparatus  we  have  employed. " He  observed  the 
resonant  enhancement  of  third  harmonic  generation  ob- 
tained by  timing  the  laser  to  the  two  photon  6S  - 9Dj^ 
resonance  and  demonstrated  that  it  was  correlated  with 


Continuum 


two  photon  absorption.  He  also  observed  intense  radia- 
tion at  two  of  the  6S  - nP  transitions  which  were  at- 
tributed to  the  parametric  frequency  conversion  process. 

This  work  should  be  considered  an  extension  of 
Leung's  observations.  We  have  made  more  detailed  mea- 
surements for  each  of  the  three  processes  emphasizing 
the  dependence  of  each  on  laser  intensity  and  laser 
frequency.  The  laser  light  source  has  been  improved  in 
bandwidth,  mode  structure,  tunability,  and  stability 
compared  with  that  used  by  Leung.  The  laser,  other 
experimental  apparatus,  and  experimental  methods  are 
described  in  chapter  2. 

Our  data  involving  third  harmonic  generation  and 
two  photon  absorption  is  presented  and  discussed  in 
chapter  4 and  the  parametric  data  and  discussion  are 
contained  in  chapter  5.  In  all  three  processes,  impor- 
tant discrepancies  are  found  between  the  observed  be- 
havior and  that  predicted  by  the  simple  theory  of  each 
which  is  presented  in  chapter  three.  Consequently,  the 
theory  is  modified  in  chapters  four  and  five  to  include 
such  effects  as  saturation  of  atomic  transitions,  photo- 
ionization, and  Stark  shifts.  We  demonstrate  that 
saturation  of  transitions  is  important  in  each  case  and 
limits  frequency  conversion  efficiency. 


CHAPTER  2 


THE  EXPERIMENTAL  APPARATUS 

The  apparatus  for  a typical  experiment  consists  of 
a tuneable  ruby  laser,  a sample  cell  containing  cesium 
vapor,  three  photomultiplier  detectors  and  their  associ- 
ated optics,  and  an  electronic  data  acquisition  system. 

The  laser  is  continuously  tuneable  in  the  range  of  694.3 
nm  to  693.5  nm.  Its  wavelength  is  determined  by  the 
temperature  of  the  ruby  and  by  etalons  in  the  laser  cavi- 
ty. Cesium  vapor  at  pressures  up  to  .5  torr  is  contained 
in  a stainless  steel  cell  fitted  with  three  windows. 

Bursts  of  light  are  detected  by  the  photomultipliers  and 
the  resulting  electrical  pulses  are  fed  to  the  data  acqui- 
sition system  which  digitizes  the  data.  A more  detailed 
description  follows.  See  fig.  2.1  for  a schematic  dia- 
gram of  the  system. 

2.1  The  Cesium  Cell 

The  cell  containing  the  cesium  vapor  is  about  20  cm 
long,  constructed  of  stainless  steel  and  fitted  with 
three  sapphire  windows. 4 It  is  lined  with  a thin  sheet  of 
copper  felt  in  the  central  portion.  These  materials  are 
chosen  for  their  resistance  to  damage  by  the  cesium.  Be- 
fore being  filled  with  cesium,  the  cell  is  evacuated  to 
about  10  ^ torr  and  approximately  one  torr  of  neon  is  ad- 
mitted to  act  as  a buffer  gas. 
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The  distribution  of  cesium  vapor  inside  the  cell  is 
determined  by  the  temperature  distribution  of  the  cell. 
Basically,  the  central  region  is  kept  hot  and  the  ends 
cool  to  produce  a peaked  distribution  about  5 cm  in  width. 
Cesium  vapor  then  circulates,  being  generated  in  the  hot 
center  and  condensing  at  the  cooler  ends.  The  condensed 
liquid  cesium  is  returned  to  the  cell  center  by  the  copper 
felt  which  acts  as  a wick.  The  neon  buffer  gas  ensures 
a gentle  circulation  and  a smooth  density  distribution 
since  its  pressure  is  higher  than  that  of  cesium. 

The  vapor  density  is  controlled  by  the  temperature  at 

the  cell  center.  The  signal  from  a thermocouple  attached 

there  is  the  input  for  a mi llivoltmeter  and  a time- 

proportioning  temperature  controller  which  maintains  cell 

temperature  within  about  1°C  of  its  set  point.  The  use- 

o 

ful  temperature  range  is  100  - 300  C corresponding  to 

-4 

a cesium  pressure  range  of  roughly  10  to  1 torr.  The 
temperature  controller  regulates  the  current  through  an 
insulated  heating  tape  wound  around  the  central  part  of 
the  cell.  The  cooling  is  provided  by  water  circulating 
through  copper  tubing  wrapped  around  the  cell  near  each 
end. 

The  cell  windows  are  made  of  sapphire,  which  is 
more  resistant  to  damage  by  hot  cesium  vapor  than  glass 
or  fused  quartz.  Sapphire  is  a uniaxial  crystal  and  its 
optical  properties  are  temperature  dependent.  To  prevent 


rf--- : ’ ~ ^ " -*  


problems,  the  end  windows  are  maintained  at  a constant 
temperature  somewhat  higher  than  the  rest  of  the  cell. 


This  ensures  that  no  condensation  on  the  windows  can 
occur  and  that  thermally  induced  changes  in  the  optical 
properties  are  minimized.  In  addition  the  windows  are 
oriented  so  that  their  optic  axes  lie  in  the  plane  of 
polarization  of  the  laser  light  so  that  the  polarization 
of  light  entering  and  exiting  the  cell  is  unaltered. 

Since  the  integral  of  the  cesium  density  along  the 
length  of  the  cell  is  the  quantity  of  interest  more  often 
than  the  density  at  the  cell  center,  the  former  is  cali- 
brated against  cell  temperatures.  This  is  done  by 
measuring  THG  versus  cell  temperature  with  the  laser 
tuned  near  the  cesium  two  photon  resonance.  A plot  of 

■ 

this  curve  is  shown  in  fig.  2.2.  As  will  be  shown  in 
chapter  3,  the  signal  at  third  harmonic  is  proportional 
to  sin 2(/  — dx'jwhere  Ak,  the  linear  dispersion  of 

V o 2 / 

cesium  between  third  harmonic  and  laser  frequencies,  is 
proportional  to  cesium  density.  Thus  the  data  of  fig. 

2.2  give  a measure  of  the  relative  integrated  density 
directly  in  coherence  lengths  where  one  coherence  length 
(1  CL)  is  defined  as  the  value  of  /Ndx  for  which 
f—  dz  = JL.  Most  of  the  data  will  be  presented  with  the 
number  of  coherence  lengths  as  the  unit  of  /Ndx.  The 
conversion  from  coherence  lengths  to  absolute  density 
can  be  found  in  two  ways.  The  first  is  to  calculate  Ak 
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and  estimate  the  density  distribution  of  cesium  in  the 
cell.  The  second  is  to  use  the  data  in  fig.  2.2  and 
estimate  the  density  and  distribution  at  a specific 
number  of  coherence  lengths.  The  two  methods  agree  well 
that 


±1 

1 CL  = .2  torr  x 5 cm  x 2 


4 x 10 
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atoms 

— r x 

cm 


5 


cm  x 


(2.1) 


2.2  The  Attenuator 

The  power  of  the  radiation  incident  upon  the  cesium 
cell  is  systematically  controlled  by  means  of  a variable 
attenuator  while  laser  output  is  maintained  at  a constant 
level.  This  ensures  that,  at  the  cell,  the  beam  charac- 
teristics such  as  time  development  and  mode  structure  are 
constant  while  the  power  is  varied.  The  attenuator  is  an 
absorption  cell,  filled  with  a solution  of  NiSC>4  in  water 
and  has  a path  length  variable  from  0-15  mm.  Trans- 
mission is  calibrated  by  measuring  power  before  and  after 
the  cell  and  is  found  to  be  T = 10~t//5’6  where  T is  trans 
mission  and  t is  absorber  thickness  in  mm. 

2.3  The  Detector  & Monitor  Boxes 

There  are  three  photomultiplier  light  detectors  and 
associated  optics  (fig.  2.1).  The  first  detector  is 
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used  to  monitor  the  shot  by  shot  performance  of  the 
laser.  For  this  a fraction  of  the  beam  is  split  off  by 
a pellicle  reflector.  The  reflected  beam  is  in  turn  split 
into  two  parts,  one  of  which  is  directed  to  the  trigger 
photodiode  to  trigger  the  data  acquisition  electronics. 

The  remainder  generates  second  harmonic  in  a quartz 
crystal.  The  harmonic  is  used  because  it  is  generally 
best  to  use  an  nth  harmonic  monitor  if  nth  harmonic 
generation  is  being  studied.  The  monitor  and  signal  then 
both  respond  to  fluctuations  in  the  laser  pulse  shape  in 
the  same  way  so  there  is  less  scatter  or  systematic  change 
in  the  ratio  of  signal  to  monitor.  Third  harmonic  genera- 


tion, two  photon  absorption,  and  parametric  gain  are  of 
interest  here  so  no  single  monitor  is  ideal.  The  second 
harmonic  one  is  convenient,  however,  so  it  is  used. 

Since  there  are  no  systematic  changes  in  pulse  shape  this 
will  not  affect  our  results  except  to  produce  some  ad- 
ditional scatter.  The  fundamental  is  removed  by  three 
centimeters  of  CuSO^  solution  and  any  other  light  is  re- 
jected by  a narrow  band  interference  filter  leaving  only 
the  second  harmonic  to  be  detected  by  the  RCA  model  1P28 
monitor  photomultiplier.  Its  signal  is  fed  to  the  moni- 
tor channel  of  the  data  acquisition  system. 

The  fundamental  or  laser  beam  is  next  attenuated  by 


the  variable  attenuator  and  in  some  cases  focused  into 
the  center  of  the  sample  cell  by  an  18  cm  focal  length 
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lens.  Light  enters  and  exits  the  ce) 1 through  the  three 
sapphire  windows  which  are  transparent  between  7000  and 
200  nin . 

The  side  detector  located  at  the  side  window  of  the 
sample  cell  is  used  to  measure  the  amount  of  fluorescence 
at  585  nm.  The  box  contains  a collecting  lens,  a CuSO^ 
filter,  a neutral  density  filter,  a narrow  bandpass 
interference  filter  at  585  nm  and  a 1P28  photomultiplier 
tube.  Figure  2.3  is  a diagram  of  this  detector  system. 

The  third  detector  looks  at  the  third  harmonic  sig- 
nal, one  of  the  strong  lines  at  the  6S  to  nP  resonances, 
or  any  other  radiation  within  the  detector  bandwidth 
(700  nm  to  200  nm) . The  wavelength  is  selected  by  tuning 
a small  quartz  prism  monochromator.  The  detector  box 
contains  an  assortment  of  filters  such  as  CuSO^,  NiSO^, 
neutral  density,  and  interference  filters  which  is  appro- 
priate for  the  signal  of  interest.  The  phototube  is  a 
1P2  8 . 

A different  third  detector  (not  shown)  is  used  when 
infrared  signals  are  under  study.  This  consists  of  a 
gold-doped  germanium  detector  operated  at  liquid  nitrogen 
temperature,  an  interference  filter  at  the  proper  wave- 
length, and  colored  glass  filters  to  eliminate  the  laser 
light . 


Fluorescence  detection  system 
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2.4  The  Data  Acquisition  System 

There  are  three  channels  like  the  one  drawn  in  fig- 
ure 2.1  which  share  the  gate  generator.5  The  pulses 
generated  by  the  four  photodetectors  are  routed  via  co- 
axial transmission  lines  to  the  electronics.  The  signals 
from  the  three  photomultipliers  are  fed  to  three  inte- 
grators and  the  fourth  signal,  the  trigger,  causes  the 
gate  generator  to  produce  a pulse  of  preset  duration  ad- 
justed to  span  the  signal  pulses.  This  pulse  gates  the 
integrators  and  strobes  the  ADCs . The  voltage  from  each 
integrator  is  proportional  to  the  charge  from  its  photo- 
multiplier and  hence  to  the  number  of  photons  for  that 
channel.  The  ADC  converts  this  from  an  analog  to  a 
digital  form  and  the  result  is  typed  out  on  a teletype. 
The  net  result  is  that  each  laser  shot  generates  a print- 
out of  three  numbers  proportional  to  the  three  signals . 
These  data  are  then  analyzed  in  a manner  appropriate  to 
the  particular  experiment. 

2.5  The  Laser 

The  laser  is  a Spacerays  Model  1010c  Q-switched  ruby 
laser  operated  at  about  1 MW  peak  power  with  pulse  dura- 
tion of  twenty  to  fifty  nanoseconds  and  a beam  area  of 
2 

.03  cm  . The  output  has  a narrow  spectral  width  and  is 
continuously  tunable  over  the  range  of  the  two  photon 
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resonance.  Figure  2.4  is  a schematic  outline  of  the 
laser  cavity. 

The  wavelength  at  which  the  laser  operates  is  de- 
termined by  three  factors:  the  temperature  of  the  ruby 
amplifying  medium,  the  optical  thickness  of  the  front 
cavity  reflector,  and  in  some  cases  by  the  tuning  of  an 
intracavity  etalon.  The  frequency  dependences  of  these 
elements  are  shown  in  fig.  2.5. 

The  temperature  of  the  ruby  rod  tunes  the  laser  be- 
cause the  fluorescence  peak  of  the  lasing  transition  in 
ruby  becomes  narrower  and  shifts  to  shorter  wavelength 
as  the  crystal  is  cooled.  Since  gain  is  large  enough  to 

support  lasing  only  near  the  peak,  the  lasing  frequency 

-1 

is  determined  within  . 5 cm  or  so  by  ruby  temperature. 

The  lasing  wavelength  near  150°K  is  given  by  A,  = 

lasing 

693.54  + .0029  (T  - 150)  where  the  unit  of  A is  nm  and 
that  of  T is  degrees  Kelvin.6 

Ruby  temperature  is  regulated  by  varying  the  rate 
of  flow  of  liquid  nitrogen  past  the  rod.  A temperature 
controller  senses  the  temperature  by  a thermocouple 
attached  to  the  rod  holder  and  opens  or  closes  a solenoid 
valve  in  the  liquid  nitrogen  line  as  needed.  In  this 

• • 

manner  the  temperature  at  the  time  of  firing  the  laser 
is  maintained  within  about  1°K  of  the  set  point  of  the 
controller.  The  temperature  actually  cycles  through  a 
range  of  3 - 5°K  with  a period  of  2 - 3 minutes  and  the 


Laser  optical  cavity  and  tuning  system 
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Fig.  2.5 


Frequency  dependence  of  laser  tuning  elements. 
Note  scale  expansion  by  10  in  upper  half. 


21 


I 


laser  is  fired  as  it  drifts  through  the  proper  tempera- 
ture. The  rod  temperature  necessary  for  resonance  with 

o 

cesium  is  about  150  K. 

The  second  tuning  element  is  the  front  reflector  of 
the  laser  cavity.  It  consists  of  a temperature  con- 
trolled, 1.0  cm  thick,  fused  silica  window  with  plane 
and  parallel  uncoated  surfaces  having  the  reflectivity 
indicated  in  fig.  2.5.  Lasing  is  found  to  occur  only 
within  a small  frequency  range  near  the  peaks  of  reflec- 
tivity and  on  only  one  such  peak  if  that  peak  is  near  the 
ruby  fluorescence  maximum.  In  this  case,  usually  5 or  6 
longitudinal  modes  of  the  laser  cavity  participate,  and 
the  intermode  spacing  is  125  Mhz . The  ruby  temperature 
is  stable  enough  that  lasing  on  only  a single  selected 
reflection  peak  is  assured. 

The  positions  of  the  reflectivity  peaks  are  tuned 
by  varying  the  temperature  of  the  reflector  and  hence 
its  optical  thickness.  Using  the  values  of  = 

.06  x 10_6/°C  and  i ^ = 7 x lo"6/°C,  where  L is  thick- 
ness  and  n is  index  of  refraction,  the  tuning  rate 
is  determined  to  be  one  order  per  2.5°C  change  in  tem- 
perature. Since  temperature  stability  is  about  .05°C 
the  stability  of  the  laser  output  is  about  .005  cm 
with  a width  of  about  .02  cm  \ The  relative  distribu- 
tion of  power  among  the  longitudinal  modes  is  not  con- 
stant, however,  so  effective  stability  is  ^ .01  cm 
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The  tuning  range  is  about  1 cm  Most  of  the  data  was 
taken  with  these  conditions. 

To  provide  continuous  tunability,  the  reflector  and 

fluorescence  peaks  are  brought  into  coincidence  then  tuned 

together  to  the  desired  frequency.  The  expected  tuning  is 

verified  by  examining  the  laser  output  with  a Fabry-Perot 

interferometer.  Figure  2.6  exhibits  the  resulting  inter- 

2 

ferograms  and  a plot  of  the  value  of  R (where  R is  the 
radius  of  the  ring  in  the  photos)  versus  the  ruby  tem- 
perature. Knowledge  of  the  plate  separation  allows  a 
determination  of  the  tuning  curve  in  cm  V °C  from  the 
graph.  The  vertical  lines  and  the  sharper  of  the  two  ring 
patterns  in  each  interferogram  are  produced  by  a helium- 
neon  laser  used  to  stabilize  the  Fabry-Perot.  It  is  clear 
from  the  sloping  lines  in  fig.  2.6  that  tuning  is  well 
behaved.  Similar  interferograms  reveal  the  spectral 
width  and  number  of  longitudinal  modes  in  the  laser  out- 
put. 

The  third  tuning  element,  the  intracavity  etalon, 
was  used  when  a narrower  spectral  bandwidth  was  desired. 
When  in  place  it  allowed  only  one  or  two  longitudinal 
modes  to  lase.  Stability  was  sufficient  to  select  any 
mode  within  the  tuning  range.  Its  insertion  into  the 
cavity  lowered  the  power  output  by  a factor  of  five  or  so, 
however,  and  the  tuning  range  was  decreased  to  about 
.05  cm  ^ . 


Tamp 


Fig.  2.6 


Top:  Fabry-Perot  inter ferograms 

Bottom:  Frequency  - ruby  temp,  calibration 


24 


This  intracavity  etalon  is  a solid,  fused  silica,  two 
surface  etalon  with  dielectric  coatings  of  80%  reflecti- 
vity and  a thickness  of  1.47  cm.  The  finesse  is  about  15. 
It  allows  lasing  only  at  its  transmission  peaks  since 
reflected  light  is  lost  from  the  laser  cavity.  Tuning  is 
achieved  by  varying  the  angle  of  tilt  of  the  etalon  or  by 
changing  its  temperature.  In  practice,  rough  tuning  is 
by  tilt  and  fine  tuning  by  temperature  control.  For  the 
latter  the  temperature  is  regulated  to  .01°C  which  is 
sufficient  to  pick  any  desired  laser  mode. 

Both  the  Fabry-Perot  and  a fast  oscilloscope- 
photodiode  combination  (Tektronix  519  oscilloscope  and 
ITT  FW  114A  photodiode  with  a combined  rise  time  less 
than  .5  ns)  were  used  to  study  the  output  when  the  etalon 
was  in  place.  The  Fabry-Perot  interference  patterns  re- 
vealed that  there  were  only  one  or  two  modes  lasing,  and 
study  with  oscilloscope  of  the  modulation  of  laser  power 
due  to  interference  between  longitudinal  laser  modes 
verified  this  and  showed  that  a single  mode  could  be  made 
to  predominate,  thus  effectively  producing  a single  mode 
laser  with  a power  of  100  kW. 

Unfortunately,  this  single  mode  operation  was  of 
limited  usefulness  because  of  the  concomitant  sacrifice 
in  power  and  tuning  range.  The  power  was  not  sufficient 
for  most  of  the  experiments  and  the  limited  range  made  it 
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CHAPTER  3 


BASIC  THEORY  OF  TWO  PHOTON  ABSORPTION, 

THIRD  HARMONIC  GENERATION  AND  PARAMETRIC  GAIN 

This  chapter  presents  the  basic  theory  for  the  three 
nonlinear  processes  of  interest,  and  applies  it  to  the 
specific  case  of  ruby  laser  radiation  interacting  with 
cesium  atoms.  The  primary  concerns  are  evaluating  the 
appropriate  x's  (see  definitions  in  table  3.1),  and 
solving  driven  wave  equations  in  which  the  induced  polari 
zations  are  the  driving  terms.  These  calculations  indi- 
cate the  size  of  the  signals  to  be  expected,  their  depen- 
dence upon  laser  intensity,  their  resonant  behavior,  and 
the  effects  of  dispersion.  Other  factors  such  as  hyper- 
fine  splittings,  laser  mode  structure,  absorption,  and 
saturation  are  necessarily  included  in  a full  description 
of  the  processes,  but  a discussion  of  these  is  postponed 
until  chapter  4. 

Before  proceeding  with  the  theory  of  the  individual 
processes,  some  discussion  of  the  characteristics  of 
the  x's  is  in  order.  It  was  mentioned  in  the  introduc- 
tory chapter  that  x is  a tensor  and  depends  on  the  fre- 
quencies of  the  applied  fields.  Taking  a third  order  X 
as  an  example,  the  polarization  is  written  in  the  form 
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to  account  for  the  tensor  nature  and  frequency  dependence 
of  X*  The  conventions  are  the  following:  is  the 

amplitude  of  polarization  in  the  direction  i(i=fc,  £,  or  2) 
at  the  frequency  w produced  by  the  three  electric  fields 
on  the  right  hand  side;  w = + K is 

listed  in  the  definitions  in  table  3.1;  the  E's  are 


Fourier  components  of  the  applied  field  at  the  frequencies 
indicated  by  the  superscript  and  with  the  polarization 
directions  given  by  the  subscripts;  each  combination  of 
Fourier  components  has  an  element  of  x associated  with  it. 
In  an  isotropic  medium  such  as  cesium  vapor,  sym- 


metry allows  only  those  elements  of  the  tensor 

which  are  of  the  form  x . • • . > X ■ ■ ■ ■ • X ■ ■ ■ ■ > or  X • 

un  ii]]  1313  ID  3i 

(j,  i=x,  y,  z;  i*j)  to  be  nonzero.  In  addition  the  fol- 


lowing relation  holds 


The  elements  of  x f°r  an  atomic  gas  can  be  derived 
from  the  atomic  susceptibility,  xat°m»  by  averaging  over 
angular  orientations  of  the  atoms,  applying  local  field 
correction  factors,  and  averaging  over  atomic  velocities 
(Doppler  averaging) . Cesium  atoms  are  symmetric  so  the 
orientational  averaging  has  no  effect,  and  the  fields  at 


u 
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Table  3.1 


Definitions 


i r -i 

= j E e 


-i(tot-kx)  -to 
+ E e 


i (ut-kx)' 


E = E 


E is  an  electric  field 


p = 1 ’pw  e-i  (ut-kx)  + p-uei  (ut-kx)" 


tO*  -GO 

p = p 


P is  a polarization 


to  to  to 

toa  12m 

P.  = Kx,.,  (-to  ; to.,  to., ..to  ) E.  E . . .E 

l i jk . . . r a 1 2 mjk  r 

K = 2nD 

n = number  of  nonzero  P^  labels  minus  number  of  non- 

zero Ew  labels 

D = number  of  distinguishable  arrangements  of  field 

to ' s . ie.  , of  to,  ...  u) 

1 m 

ij...r  refer  to  polarization  directions  St,  and  2 


Table  3.1  (cont.) 


I 

Complex  numbers  or  functions  are  written  as:  x = x + 
ix"  where  x'  and  x"  are  real. 


I • - •• 


Z(x,y)  is  the  plasma  dispersion  function 

00  2 

Z (x,y)  = — fm  { e_t  /(t  - X - iy)  } dt 

/tt 


‘n£ 


CO 


n£ 


_i  Fn& 


CO 


n£ 


is  the  frequency  of  the  transition  from  state  n 

to  state  SL  and  F „ is  the  FWHM  for  that  transition. 
n£  . 

If  state  £ is  the  ground  state,  F2  0 is  written  fi  . 

nx,  n 


03 


0)  = 

2ttc 

A<2n£ 

= to. 

power  of  radiation  at  <o 
Q = 

2\h 


<o 


n 

2ttc 


n£ 


{2kT/Mc}’ 


Doppler  width 
2/£n2 


the  atomic  sites  are  just  the  applied  fields  if  the  gas 

is  dilute  as  it  is  in  our  case.  Therefore,  y is 

AijkX, 

simply  N • <x . . > where  N is  the  number  of  atoms 

ljkJl  Doppler 

per  unit  volume. 

The  susceptibilities  which  appear  in  the  remainder 

of  this  text  are  yatom  unless  indicated  otherwise  so  the 

zzzz 

superscript  and  subscripts  will  be  omitted  for  brevity. 
Doppler  averaging  will  be  indicated  by  <x>* 


3.1  Two  Photon  Absorption  (TPA) 


In  this  section  an  expression  relating  the  nonlinear 
susceptibility  x to  a two-photon-absorption  cross  section 
is  derived  by  analogy  with  the  one  photon  absorption  case. 
The  numerical  value  of  the  cross  section  is  given,  Doppler 
averaging  is  performed,  and  the  hyperfine  splitting  of  the 
ground  state  is  discussed. 


Just  as  linear,  or  single  photon,  absorption  is  pro- 
duced by  the  imaginary  part  of  the  linear  refractive 
index,  two  photon  absorption  is  caused  by  an  imaginary 
part  of  the  nonlinear  refractive  index.  Specifically,  it 
will  be  shown  below  that  two  photon  absorption  depends  on 
X"  (-0) ; oj,  u), -w)  while  single  photon  absorption  depends  on 
X"(-w;w)  where  x"  is  the  imaginary  part  of  x (see  table 
3.1)  . 

If  all  contributions  to  the  polarization  at  fre- 
quency u),  created  by  a single,  linearly  polarized  applied 
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field  of  frequency  u>  are  considered,  the  polarization  in 
an  isotropic  medium  may  be  written  as 


CO  (0  (0  (0  — (0 

P = NK^X-^-W/W)  E + NX^Xj  _w)  E E E 

(3.3) 

+ NK5X5  (“<*>;  <*> , u),  -co,  -a),  co)  EW  E W E W E^  EW  + ... 


Keeping  only  the  first  two  terms  we  have 


Pw  = n{KiXi  ( — u) ; (jo)  + K3x3(-w;u),to,-oj)  Ew  E-liJ}  Ew 

(3.4) 

The  first  term  accounts  for  the  linear  index  of  refrac- 
tion and  single  photon  absorption.  It  is  useful  to  con- 
sider the  second  term  an  intensity  dependent  extension  of 
the  first.  This  will  produce  an  intensity  dependent  re- 
fractive index  and  an  intensity  dependent  absorption  cross 
section.  The  term  inside  the  brackets  is  now  to  be 
treated  as  the  effective  susceptibility.  Substituting  it 
in  the  usual  definitions  of  the  dielectric  constant  (e) 
and  the  refractive  index  (n) , we  have 


£ = 1 + 4ttN 


r CO  “(0 1 

{ K1X1  + k3x3  E e } 


(3.5) 


. — ( CO  “C0 1 

n = /e  = 1 + 2ttN  {KlXl  + K3x3  E E } 


(3.6) 
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In  this  last  step  it  was  assumed  (e-l)<<l.  it  is  now 
possible  to  find  the  absorption  cross  section. 

2 

o = -y-  Im{K1x1  + K3x3  Ew  E_tJ } = o;L  + a (3.7) 

2 2 
8 TT  m 8tT  ••  (jj  ~co 

°1  = X K1X1  °2  = A K3x3  E E (3.8) 

The  cross  section  cy,  which  produces  linear  absorption, 

II 

is  independent  of  intensity  and  proportional  to  x1(-(JJ;aJ)  • 
The  cross  section  c^,  on  the  other  hand,  depends  on 

It 

X3  (-<jo;o),  id,  -w)  and  on  the  laser  intensity,  so  the  proba- 
bility that  an  atom  will  absorb  a photon  due  to  o3  is 
proportional  to  the  intensity  squared,  as  is  expected  for 
TPA. 

II 

In  order  to  find  the  value  of  o^,  X3  must  be  evalu- 
ated for  a cesium  atom  in  the  ground  state  (6S,  ) . Keep- 
ing  only  the  term  which  is  nearly  two  photon  resonant,  we 
find  that  when  all  fields  are  polarized  in  the  z direc- 
tion2 


X3  (-w;  -w)  = (-h)  3 2R1  2 3 e4 


(3.9) 


<6S,  | z|  fcP><JlP  | z|  9D,  ><9D,  | z|  nP><nP|  z | 6S  > 
)'  2 H 1 

nf,  ic 

(fiip-  0))  (n  - u)  (n  - 2u) 


* 
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where  R denotes  the  average  obtained  by  permuting 


the  E fields  and  ft,,  is  ft~^ 

9 9D3-6S, 

this  simplifies  to 


x3  - e4  <-W'3  Z (fi„p-  2“rl 


Since  nlp  - ntp 


(3.10) 


i <9Ds I z I nP><nP I z I 6S>  > 2 

I — 5 — r : - 


(ft  -Cl)) 
nP 


This  has  been  evaluated  by  Orr  as 


X3  * 2.25  x 10-31  (ft  - 2 w)"1  esu 


(3.11) 


where  the  units  of  ft  and  to  are  cm  . This  is  the  value 
of  X3  for  a single  atom  of  cesium  in  the  6S^  state.  For 
comparison  with  experiments,  Doppler  averaging  must  be 
performed.  This  yields 


= 2.25  x 10 


-31  Z ( x, y ) 


(3.12) 


where  Z(x,y)  is  the  plasma  dispersion  function,9  Adig  is 
the  Doppler  width  in  cm  3 divided  by  2/TnJ,  and 


„ • _ " 

2d5-ftq  -ftq  ^ , 2 it 

x = ; y = — - ; A(ig  = ftg{2kT/Mc  P (3.13) 

Auig  Au>9 
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Taking  the  imaginary  part  of  <X3>  and  using  Aoig  = 

.025  cm  gives 

<X3>  = 9.3  x 10-30  Z " (x,y)  (3.14) 

Finally,  using  equations  3.8  and  3.16,  we  find 

o 2 = 6.7  x io“26  iz" (x,y)  (3.15) 

2 o 

where  I is  in  watts/cm  and  a is  in  cm^ . It  should  be 

II 

mentioned  that  the  maximum  value  of  Z is  approximately 

/if  and  occurs  when  the  laser  is  tuned  to  half  the  6S.  -9D3 

H -s 

transition  frequency  (x  = 0) . The  resonance  shape  of  a 2 

II 

is  shown  in  fig.  3.1  where  Z is  plotted  as  a function  of 
x for  y = .1.  Note  that  according  to  its 
definition  in  eqn.  3.13,  y is  approximately 
the  ratio  of  the  homogeneous  width  to  the  Doppler 
width. 

The  expression  for  03  is  modified  if  the  hyperfine 
splitting  of  the  ground  state  (6S^)  is  included.  The 
state  is  split  into  two  hyperfine  components  with  a 
separation  of  .31  cm-'*"  and  the  statistical  weights  7/16 
and  9/16  so  eqn.  3.15  becomes 


0.1 


0.1  em"1 


Fig.  3.1 

Top:  Z"  (x, . 1) 

Bottom:  hyperfine  split  TPA  cross  section 

\ 
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f,  = (7/16)  6.7  x 10-26  iz"{x  - h,  y } + 

^ 16 


16 


(3.16) 


(9/16)  6.7  x 10-26  I Z " { x + _Z  h,  y} 

16 


where  h = 


.31 


ACOr 


In  summary,  the  rate  at  which  cesium  atoms  are  pro- 


moted from  the  ground  state  to  the  9D3  level  is  propor- 

z 


tional  to  the  square  of  the  laser  intensity  and  exhibits 
a resonant  behavior  described  by  eqn . 3.16  where  x is  the 
normalized  detuning  of  the  laser  from  the  two  photon 


resonance  between  the  states  65^  and  9D^,  and  y (see  eqn, 


3.13)  is  much  less  than  unity 


3.2  Third  Harmonic  Generation  (THG) 


This  section  presents  some  aspects  of  the  theory  of 


two-photon-resonant  third  harmonic  generation  in  cesium 


2 


vapor.  In  particular,  we  will  derive  the  sin  (7rNCL/2) 


dependence  of  the  third  harmonic  signal  which  we  used  in 
chapter  2 to  calibrate  the  amount  of  cesium  vapor  in  the 
cell  (see  fig.  2.2);  we  will  discuss  the  resonance  line 
shape,  including  Doppler  averaging,  since  it  will  be 
important  in  the  next  chapter.  The  presentation  of 


theory  in  this  section  follows  that  of  Leung. 


. _ 2 


Dispersion  in  the  cesium  gas  causes  the  sin  ( ttN  /2) 

CLj 


dependence  of  the  third  harmonic  signal.  Here  N , the 

VwIj 


number  of  coherence  lengths  of  Cs,  is  proportional  to 
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the  product  of  the  integral  of  the  cesium  density  along 
the  length  of  the  cell  and  AkQ,  which  is  a measure  of 
the  dispersion.  In  order  to  derive  this  result  we  con- 
sider first  a sample  of  uniform  density  N and  length  L 
in  which  we  wish  to  generate  third  harmonic.  The  third 
harmonic  polarization  is  locked  in  space  and  time  to  the 

(l) 

applied  field  E so  it  has  the  wave  vector  3k  , where  k 
is  the  wave  vector  for  Ew  in  the  sample.  The  third 
harmonic  radiation,  in  contrast,  has  the  wave  vector 
k^  in  the  sample.  Since  all  the  k vectors  point  in  the 
x direction,  we  drop  the  vector  notation  and  define  Ak  as 
3k  - k_  . A nonzero  value  of  Ak  means  that  the  phase 
velocities  of  laser  radiation  and  its  harmonic  are  dif- 
ferent in  which  case  the  third  harmonic  radiation  genera- 
ted at  one  point  in  the  sample  will  be  out  of  phase  with 

that  produced  at  a point  farther  along.  This  interference 

2 

in  the  harmonic  waves  causes  the  sin  dependence  of  THG . 
This  result  is  derived  below  by  solving  the  driven  wave 
equation . 

The  third  harmonic  polarization  P3u)  is  given  by 

3 i 3k  x 

P W(x)  = 1/4(EW)  N<x  (~3w;  a),  w,  uj)  > e (3.17) 

The  wave  equation  for  the  third  harmonic  is 

VxVxE3w(x)  - k^  E3w(x)  = 4tt  ( 3oo/c)  2 P3W(x) 

3U) 


(3.18) 
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where  the  polarization  P w(x)  acts  as  the  driving  term. 

ik^x 

If  a solution  of  the  form  E3uJ(x)  = A(x)e  w is 
assumed,  then  substitution  into  eqn . 3.18  gives 


ik,  x 


.2  **12-  (ik,  ) e 
3x  3w 


ttN  (3co/c)  ^ (EW) '^<x>  e 


(3.19) 


i3kwx 


where  3 ^.jx^  ^as  been  neglected  because  it  is  assumed  to 
3 X l; \ I 


be  small  compared  with  2k ^ 
This  simplifies  to 


3A  (x) 


iAkcNx 


(3.20) 


u 3 

9A(x)  2 (E  ) x 

— = ttN(3w/c)  -r-r <X>  e 

3x  1K3w 


where  AkQ  = -jy. 

Now  assuming  A(0)  = 0 and  integrating  from  x = 0 to  x = L 


yields 


iAkQNx 

A (L)  = C(EV  (e  Ak  N ~~  <X>  N 

o 


(3.21) 


The  power  at  third  harmonic  is  proportional  to  |A(L) 
since  A is  a field  strength,  so 


r • 
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E3w  - Sin2(AkoNL/2)  { | EW  | 2 } 3 |^|2  (3.22) 

This  is  the  sin  dependence  mentioned  earlier. 

The  density  of  cesium  was  assumed  uniform  in  deriv- 
ing eqn.  3.22.  This  is  not  a good  approximation  to  the 
experimental  case  but  eqn.  3.22  is  easily  extended  to  the 

case  of  nonuniform  density  by  replacing  the  arguement  of 
Ako  L 

the  sine  by  — — fQ  N(x)dx  where  N (x)  is  the  cesium  number 
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It  is  interesting  to  note  that  both  the  real  and 
imaginary  parts  of  x contribute  to  the  generation  of 
third  harmonic.  The  harmonic  generated  by  the  imaginary 
part  of  x is  90°  out  of  phase  with  that  produced  by  the 

,2  m2  ,2  .,2 

real  part  so  they  add  as  x + X or  Z + Z Par 

I 

from  the  resonance  THG  is  due  solely  to  Z while  at  the 

ll 

line  center  only  Z contributes.  The  total  lineshape 
2 2 

Z + Z has  the  resonance  shape  shown  in  fig.  3.2  with 

I 

the  resonance  peak  at  2m  = or  x = 0. 

If  the  hyperfine  splitting  of  the  ground  state  is 
to  be  included,  the  Doppler  averaged  x will  have  the  form 


<X  (-3uj;  oj,  gj,gj)  > « 


(3.25) 


( 7/16 ) Z { x - yl  h,  y}  + ( 9/16 ) Z { x + h,  y} 


• 31  2 

where  h = Since  is  proportional  to  |x|  it  must 

have  the  form 


(3.26) 


2 

[(7/16)z'{x  - x,  y}  + (9/16) Z ' { x + h,  y}]  + 


[(7/16)z"{x  " xf  x,  y}  + (9/16) z" { x + jL  h,  y}] 


This  is  graphed  in  fig.  3.2. 


0.1 


- 0.1  cm' 


Fig.  3.2 

Top:  Frequency  dependence  of  Z'(x,.l),  Z" (x, . 1) , | Z (x, . 1) | 

Bottom:  hyperfine  split  THG  efficiency 
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when  the  laser  is  tuned  between  the  two  resonances 

I 

the  Z terms  have  opposite  signs  whereas  they  have  the 
same  sign  if  the  laser  is  higher  or  lower  than  both.  This 
interference  produces  asymmetry  in  the  two  resonance 
peaks  which  moves  their  centers  of  weight  apart  somewhat. 

II 

The  Z terms  are  narrow  in  relation  to  their  separation, 
and  their  interference  does  not  significantly  affect  the 
lineshape . 

To  summarize,  the  characteristics  of  the  third  har- 
monic process,  as  predicted  by  the  theory  discussed  above, 

2 

are:  linear  dispersion  produces  a sin  dependence  on  the 

total  amount  of  cesium  traversed;  the  signal  varies  as  the 
cube  of  the  laser  power  if  the  number  of  coherence 
lengths  of  cesium  and  laser  tuning  are  constant;  the 
resonance  shape  is  the  Doppler  broadened  resonance  line 
shape  as  shown  in  fig.  3.2. 

3.3  Parametric  Gain 

In  this  section  we  present  the  theory  for  the  reso- 
nantly enhanced  parametric  process.  The  section  is 
divided  into  several  subsections,  the  first  of  which  con- 
tains a brief  derivation  of  the  equations  describing 
parametric  gains.  The  resonant  behavior  of  gain  is  de- 
duced from  these  equations  in  the  second  subsection;  and, 
in  succeeding  subsections,  the  existence  of  a threshold 

power  and  numerical  estimates  of  gains  and  threshold 
powers  are  discussed. 


_i 
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3.3a  The  Gain  Equations 

The  parametric  process  may  be  considered  to  be  two 
separate  three-wave-mixing  processes  which  produce  re- 
spectively the  signal  and  idler  polarizations.  The  first 
of  these  mixes  the  pump  (laser)  wave  with  itself  and  with 
the  idler  wave  via  x (-ws ' ' WP ' -wi > * T^e  second  mixes  the 

signal  and  pump  waves  via  x (~wT ; w # w , -u) „)  . The  polariza- 

JL  ir  ir  O 

tions  at  the  signal  and  idler  frequencies  are  in  turn  the 
source  of  radiations  which,  depending  on  their  phases 
relative  to  the  signal  and  idler  waves,  can  produce  gain 
or  loss  in  the  signal  and  idler. 

The  derivation  of  the  gain  equations  presented  in 
this  subsection  follows  that  of  Yariv  and  Pearson.10  The 
application  to  the  particular  case  of  the  ruby  laser- 
cesium  vapor  has  been  developed  by  Leung.4  The  deriva- 
tion is  displayed  here  in  order  that  approximations  may 
be  pointed  out  and  to  facilitate  the  discussion  of  reso- 
nance lineshapes  in  a later  subsection. 

From  the  description  of  the  parametric  process,  it 
is  apparent  that  in  order  to  find  the  intensity  of  the 
idler  and  signal  fields,  we  must  solve  a pair  of  coupled 
wave  equations  in  which  the  polarizations  at  the  signal 
and  idler  frequencies  act  as  the  source  or  driving  terms . 
The  polarizations  are  given  by 
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i(kD  + k„  - k)  x 


S p p _T  T ~ X,. 

P (x)  = NKg<xs>E  E E 1 e P P 1 


PMx)  = NKi<Xi>Ef  Ef  E-s  e 


P pP  T?-S  ^ ^kp  + kp  ks^X 


(3.27) 


where  xs  = x(-“s.*up#up#-uI)  and  Xj  = X (-Wj ; wp , wp,  -u  ) . 
Note  that  k is  not  confined  to  real  values  so  linear 
absorption  can  be  included. 

The  coupled  wave  equations  which  are  to  be  solved 


o c -i^ct  - k x)  2 -i  (a)  t - k x) 

V (E  e ‘ ) - ^ 7T  {ES  e 

c 9t 


p p -i  “i(V  " 2kpx  + kix) 

4ttNKs<xs>E  E E e } 


(3.28) 


2 -I  i(wIt  ' kIx)  I t i(u)rt  " krx> 

V2(E  1 e ) = h H 1E  e + 

c 9 tr 

* -p  -p  s i(wxt  - 2kpx  + kgX) 
4ttNKi<xi>E  E E e } 


We  have  assumed,  in  writing  this  pair  of  equations  that 
kp  is  real  and  also  that  there  is  little  loss  from  the 
pump  wave  due  to  conversion  of  its  energy  into  signal  and 
idler  waves . 

If  the  field  strengths  EP  and  EE  are  assumed  to 
change  by  a very  small  percentage  per  wavelength,  then  the 


vfl 
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terms  involving  — j are  much  smaller  than  those  involving 
_ and  may  be  neglected.  We  are  then  left  with  the  fol- 
lowing pair  of  coupled  first  order  differential  equations. 


3x 


i2TTU) 

§-  <Xo>N(3/4) 

k„c  S 


iAkx 

e 


-i27TU>, 


v2 


<Xj>N ( 3/4) 


-iAkx 

e 


(3.29) 


where  Ak  = 2kp  - kg  - k*  . 

Equation  3.29  may  be  solved  by  conventional  methods 
to  yield  the  following  solutions  for  ES  and  E1  . 


E 


-I 


where 


{ Eq  (Coshyx  - Sinhyx)  + EoI  Sinhyx}  eislAkx 

{EoI(Coshyx  + Sinhyx)  - — y Sinhy z } e-J5lAkx 

(3.30) 


C p p 

A = — j — N<X  >E 

c 4 k S 

S 

2tt  3 a)  2 

B = -y N<x*>E“P  E-P  (3.31) 

c 4 k 1 

I 

G = AB  Y = {G  - (Ak/2)  2 } ** 
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9 

These  are  the  basic  equations  describing  the  para- 
metric process.  To  arrive  at  this  result  we  have  assumed 
there  is  little  loss  of  the  pump  intensity  due  to  either 
linear  absorption  or  conversion  to  signal  and  idler  waves 
although  linear  absorption  of  signal  and  idler  radiation 
is  allowed.  We  have  also  assumed  that  the  gain  per  wave- 
length is  small  and  have  included  Doppler  averaging  in 
eqn . 3.27  as  we  must,  rather  than  in  eqn . 3.30  where 
the  product  XgXj  would  appear. 

3.3b  Resonant  Shape  of  Parametric  Gain 

This  subsection  is  devoted  to  a discussion  of  the 

resonant  behavior  of  parametric  gain.  Since  this  is  a 

complex  matter,  we  shall  consider  first  the  nature  of  the 
* 

product  <XS><XI>  which  occurs  in  y (eqn.  3.31),  then  the 
behavior  of  Ak  will  be  discussed,  and  finally,  the 
resonance  shape  of  the  gain  itself  is  considered. 

In  the  case  of  cesium  atoms  interacting  with  ruby 
laser  radiation,  the  value  of  each  x is  maximized  when 
the  laser  is  tuned  to  the  two-photon  resonance  and  the 
signal  is  resonant  with  the  6S,  - nP  transition,  where 

the  nP  state  is  any  P state  with  energy  less  than  that 
of  9D3 . This  is  evident  from  the  following  expressions 
for  <xs>  and  <XI>- 
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<Xq>  = r-  <6S,  | z|nP><nP|  z 9D,> 

b 3 (he) ^ * * 


<9D^|  z | mPXmP  | z|  6Sj5> 

m fi  „ - uL 

mP  P 


<(6g  - 2cog)  1 (finP  - us)  1> 


(3.32) 


where  nP  is  the  state  near  resonance  with  the  signal 

frequency.  Tuning  the  laser  to  resonance  maximizes 

{fig  - 2£>p}  while  the  factor  {finp  - d>s } ^ is  maximum  if 

the  signal  is  near  the  65^  - nP  resonance.  The  same  is 
_ ★ 

true  for  <Xj>  since 


3 (he) 


<6S^|z|nP>  <nP|z|9D3>* 


<9D^  | z | mPXmP  | z | 6S^> 


(3.33) 


fi  „ - uL 

mP  P 


<(fig  - 2(2p ) 1 (finp  - dig)  X> 


Performing  the  Doppler  averaging  indicated  in  eqn . 3.30 
and  3.31  yields 


< (fi  - 2w  ) 1 (fi  - <2  ) 1>  = 

9 P nP  S 


(3.34) 


Adi  • Aw„ 

nP  9 


ZU  ) - z(e  ) 
nP  9 


^nP  - S 
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and 


<(«9  " 2<V  1 ^nP  ‘ 1>  = 


zu„p>  - Z‘(V 


(3.35) 


Aw  • Aw„ 
nP  9 


CnP  - S 


where 


C = (x  + iy  ) ; x = 

9 9 9 9 


2Wp  - fig 

Aw„ 


— fi. 


; y„  = 


9 Aw, 


£ = (x  + iy  ) ; x 

nP  nP  "nP  nP 


as  - V 


-fi 


; y 


nP 


nP 


Aw  „ *“  Aw  „ 

nP  nP 


(3.36) 

and  Z (£)  is  the  Plasma  Dispersion  Function.9 

Combining  these  results  of  Doppler  averaging  and 
eqns . 3.30  and  3.31  we  arrive  at  the  following  expression 
for  <XS><X • 

<XS>  <XI>  = 5.06  x 1018  (AwrP  • AWg)~2  • 


<6S^  | z | nPxnP  | z | 9D,  > | • 


< 9 D^.  | z | mP><mP  | z | 6S^> 


m 


fi  _ - w_ 

mP  P 


Z(Cnp)  - zu9> 


Z<«nP>  - 2 U9; 


^nP  ~ S 


4nP  ~ S 


(3.37) 
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★ 

The  resonance  shape  of  <x  > <XT>  is  just  that  of  the 
last  two  terms  of  eqn.  3.37.  Graphs  of  the  real  and 
imaginary  parts  of  each  term  as  a function  of  xnp  are 
shown  in  figures  3.3  and  3.5  for  the  laser  tuned  to 
resonance  (Xg=0)  and  in  figures  3.4  and  3.6  for  the  laser 
detuned  slightly  (x^=l) . Values  of  .1  and  .05  were  cho- 
sen for  ynp  and  yg  respectively  to  facilitate  demonstrat- 
ing features  of  the  functions.  In  reality,  somewhat 
smaller  values  should  be  used  so  the  natural  linewidth 
features  in  figures  3.3  and  3.4  should  be  narrower  and 
greater  in  amplitude,  but  the  use  of  .1  and  .05  will  not 
affect  the  conclusions  we  shall  reach  concerning  the 
maximum  value  of  gain  itself. 

Figure  3.5  shows  the  real  and  imaginary  parts  of  the 
resonance  lineshapes  of  <Xg>  when  x^  = 0 . They  have 
maximum  values  of  about  two  and  widths  of  a Doppler  width. 
The  lineshape  shown  in  fig.  3.6  (Xg=l)  is  somewhat  dif- 
ferent but  the  widths  and  amplitudes  are  about  the  same. 
The  values  chosen  for  the  y 's  makes  little  difference  in 

either  figure  as  long  as  they  are  much  less  than  1. 

* 

The  plots  of  <xI> , shown  in  figures  3.3  and  3.4, 
have  the  interesting  feature  of  a tall  spike  with  the 
homogeneous  linewidth  and  a height  which  is  sensitive  to 
the  values  used  for  the  y's.  In  fact,  it  may  be  seen 
that  for  small  values  of  y^  and  y^,  the  shape  near 
xnp  = Xg  is  determined  largely  by  the  denominator 
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(x  „ - x „)  + i(y  „ + y ) since  the  numerator  is  a slowly 
nP  9 nP  9 1 

varying  function  of  x^p  compared  to  the  denominator. 

Thus  the  width  of  the  narrow  peak  is  just 


<nP  = ^nP  + ^9 


(3.38) 


or  substituting  the  definitions  of  x and  y the  width  doo 


h nP  " 

dco  = ft  + ft„ 

nP  nP  ...  9 

w9 


(3.39) 


Furthermore,  the  maximum  amplitude  is  2 /if/ (ynP+  Yg)  if  Xg 

is  near  zero,  since  the  numerator  is  about  2i/ir  for  x^ 

and  x „ near  zero.  Notice  also  that  when  the  laser  is 
nP 

detuned  by  6c0g,  the  narrow  peak  occurs  at  a detuning 

6u>  = 6co_  (u)  /o)„ ) . 

nP  9 v nP7  9 

This  line  narrowing  effect  is  reminiscent  of  the 
laser-induced  line-narrowing  effects  discussed  by  Feld 
and  Javan11  and  several  other  authors.  In  fact  the  Dop- 
pler averaging  integrals  which  occur  above  are  identical 
to  those  appearing  in  the  Feld  and  Javan  paper.  The 
laser-induced  line-narrowing  occurs  in  a system  of  three 
energy  levels  when  lasing  is  active  between  two  of  the 
levels.  If  the  small  signal  gain  from  one  of  the  lasing 
levels  to  a third  level  is  measured,  it  is  found  that 
there  are  two  dips  or  peaks  with  characteristic 


. L • • - 
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homogeneous  width  superimposed  on  the  much  broader 
Doppler  profile.  One  of  the  bumps  is  due  to  two  succes- 
sive single-photon  transitions  where  one  transition 
populates  the  common  energy  level  with  atoms  of  a se- 
lected axial  velocity  and  the  second  transition  involves 
only  those  atoms.  The  other  bump  is  due  to  two-photon 
transitions.  It  does  not  involve  velocity  selective 
population  of  a level  but  arises  from  the  nature  of 
Doppler  averaging.  An  intuitive  understanding  of  the 
effect  is  difficult  to  achieve.  The  line  narrowing  in 
our  case  is  of  the  two-photon  type. 

The  resonance  profile  of  the  product  <x>  <xt>  i-s 

b 1 

plotted  in  fig.  3.7  for  = 0 and  fig.  3.8  for  x^  = 1. 
The  narrow  peak  is  still  a prominent  feature  in  both 

cases  but  the  real  parts  are  negative  near  x = x„ . As 

n 9 

we  shall  see  in  a moment,  the  negative  sign  produces  a 
dip  in  the  gain  profile  rather  than  an  enhancement  so 
the  values  chosen  for  the  y's  only  changes  the  width  of 
that  dip  and  have  little  effect  on  the  maximum  gain  as 
long  as  they  are  much  less  than  one. 

We  consider  next  the  resonant  behavior  of  the  other 
terms  in  the  expression  for  gain.  From  equations  3.30 
and  3.31,  the  gain  coefficient  at  high  pump  powers  is 
given  by 
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gaing  (high  pump  power)  = 

2 { G - (Ak/2)2}1*-  {k " - k"} 
1 SI 


(3.40) 


It  II 

so  we  must  next  find  the  shapes  of  Ak  and  k and  k . If 

O J. 

we  assume  that  all  the  atoms  remain  in  the  ground  state, 
then  kp  and  kj  are  real  and  constant  over  the  range  of 
interest,  but  kg  exhibits  the  usual  Doppler  broadened 
resonance  behavior.  The  form  of  k is  given  by 


ko  = 3 Z ( x , y ) 

S nP  nP 


(3.41) 


where  3 is  a constant. 

We  can  now  discuss  the  resonance  shape  of  the  gain. 

Consider  first  the  case  of  high  pump  power.  Since  G in 

2 

eqn.  3.40  is  proportional  to  , the  gain  will  be  sim- 
ply 2/G  at  sufficiently  high  pump  powers  and  the  shape 
will  be  that  labeled  C = 0 in  figures  3.9  and  3.10.  The 
dip  at  the  line  center  in  fig.  3.9,  which  was  mentioned 
earlier,  has  a width  of  approximately  a homogeneous  width 
and  occurs  because  G is  negative  and  real  at  that  point 
(see  f’g.  3.3  and  eqn.  3.31). 

As  the  laser  or  pump  power  is  lowered,  the  other 

•I  II 

terms  (Ak,  k , k ) increase  in  relative  importance  and 

b J- 

the  gain  profile  assumes  the  shapes  labeled  C = 1 through 
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pump  powers.  Note  that  positive  values  represent  gain 
in  these  plots,  and  negative  values  represent  loss. 


3.3c  Existence  of  Threshold  Pump  Power 


If  the  pump  power  is  lowered  sufficiently,  the  gain 
curve  is  negative  at  all  frequencies  as  it  is,  for  in- 
stance, in  the  curve  labeled  C = 4 in  fig.  3.9.  This 
suggests  that  there  is  a threshold  pump  power  below 
which  gain  does  not  appear.  In  fact,  at  low  values  of 
Ep,  the  expression  for  gain  which  we  have  used  in  gen- 
erating fig.  3.9  (eqn.  3.40)  is  not  valid  since  the  Sinh 
terms  of  eqn.  3.28  become  comparable  with  the  Cosh  terms. 
However,  eqn.  3.28  reduces  to  equations  describing  simple 
linear  loss  at  E = 0 so  there  must  indeed  be  a threshold 

U) 

pump  power.  Furthermore,  the  threshold  power  must  have 
a value  such  that  /G  ^ Ak. 


3.3d  Numerical  Estimates  of  Parametric  Gain 


It  is  now  possible  to  calculate  the  parametric  gain 
coefficients  for  each  of  the  five  possible  resonantly 
enhanced  signal  waves  corresponding  to  the  five  inter- 
mediate P states.  First,  we  calculate  the  gain  for  high 
pump  powers  where,  according  to  eqn.  3.30,  the  gain 
coefficient  is  2 /G.  Using  the  definitions  given  by  eqn. 
3.31  we  arrive  at  the  following  expression 
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gaing  (high  pump  power)  = 


, c w to  o T 

-15  SI  2 2 * 

2.46  x 10  (Ip/A)  N <xs>  <XI> 


(3.42) 


ns  ni 


Using 


<9Ds | z I mP><mP [ z I 6SV>  _2i 

l — = 6.5  x io  esu 


m 


(f2  - oj  ) 

mP  S 


(3.43) 


in  the  expression  for  <XC>  <XT>  we  find 

O J. 


<XS>  <X*>  = Q2  F(Cnp,  Cg)(3.70  X 10  31) 


~ 2 ~ 

cj„ 

S 9 


(3.44) 


Putting  this  in  eqn . 3.42  yields 


gain^  (high  pump  power)  = 


oj  1 

{2.74  x lo“2°  — Q2  N2  ff2  F(C  , C9)} 

ws  ns  ni 


(3.45) 

h 


where  .03  has  been  used  as  the  value  for  A,  the  beam 


area,  and  S>  is  in  megawatts.  The  factor  F(£  , t,  ) is 

p nP  9 


the  product  of  the  last  two  terms  of  eqn.  3.37  (see  plots 
in  figures  3.7  and  3.8)  and  Q is  defined  by 


Q 


<6S| z|nP><nP| z| 9D,> 


(3.46) 
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E. 


f ] ^2  • 

The  maximum  value  of  Re  {F(£  , Cn)l  is  atl°ut 

r.P  9 

1.25  (see  C = 0 curve  of  fig.  3.9)  so 


gaing  (high  pump  power)  = 


13  P Q2' 

.“s 


2.07  x 10  iJ  — Q NE 


(3.47) 


9 — 4 

For  the  lOP^  line,  • (tOj/tOg)  is  3.20  * 10  and  the 
estimated  value  of  Nz  at  one  coherence  length  is 
1.75  x 10^6  (eqn.  2.1) . Finally,  expressing  Nz  in  co- 
herence lengths,  (N  ) equation  3.47  becomes 

CL 


exp {50. 8 ffpNCL} 


(3.48) 


According  to  equation  3.47  the  gains  for  the  other 

lines  can  be  found  by  evaluating  the  appropriate 

{q2  (to  /to  )}^.  These  values  are  listed  in  table  3.2  and 
X s 

the  gains  relative  to  the  lOP^  line  are  listed  in  table 
3.3.  Thus  the  equation  describing  the  behavior  of  all 
the  lines  is 


Ec  = E°  exp {50.8  *pNclR} 


(3.49) 


where  R is  the  relative  gain  from  table  3.3. 
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Table  3.2 


Line 

(^(cm 

Ugfcm  P) 

Q 

ws 

Q2 

6E\ 

17651 

11178 

. 741 

8.68 

X 

10'1 

7063 

21766 

.0930 

2 . 80 

X 

10~3 

3120 

25709 

.0644 

5.02 

X 

10-4 

9P. 

1192 

27637 

.0834 

3.00 

X 

IQ"4 

*5 

10P. 

102 

28727 

-.300 

3.20 

X 

10~4 

h 

6 P 3 

17096 

11732 

.156 

3.55 

X 

10~2 

7P3 

6 882 

21947 

.0361 

4.09 

X 

IQ'4 

8P1 

3037 

25792 

.0331 

1.29 

X 

10-4 

9pi2 

1147 

27682 

.0569 

1.34 

X 

10-4 

10P, 

75 

28754 

-.198 

1.02 

X 

10~4 

'i 

Line 

Relative 

Table  3. 

Gain  = R 

3 

W = 

IP 

(Mw) 

T 

Threshold 

6P 

92.2 

41.5 

7P 

5.24 

2.7 

*5 

00 

’V 

►r 

2.22 

0.62 

9P 

h 

1.71 

0.17 

• 

10P 

1.77 

0.06 

6P* 

18.7 

415 

OO  -^j 

>0  *X3 

2.00 

44 

1.12 

13 

9p4 

•i  1 

4.4 

1 

10P 

\ 

1 . Uu 

2.4 

n 

From  the  values  of  R it  is  evident  that  for  large 
values  of  the  pump  power  the  gain  for  the  J = 1/2  fine 
structure  component  of  each  line  is  larger  than  for  the 
J = 3/2  component,  and  that  the  gain  increases  for  de- 
creasing values  of  the  principal  quantum  number  n. 

It  should  be  mentioned  that  in  this  calculation  we 

have  considered  only  susceptibilities  of  the  form  Xzzzz 

or,  in  other  words,  we  have  assumed  all  the  fields  are 

polarized  in  one  direction.  It  is  also  possible  to 

have  the  signal  and  idler  waves  polarized  perpendicular 

to  the  pump  since  this  would  involve  x 1 s which  are 
c r Ayyzz 


allowed  by  symmetry  to  be  nonzero  (see  eqn . 3.2) 


The 


X 's  are  smaller  than  the  for  both  the  P, 

yyzz  zzzz  ^ 

and  cases,  however,  so  the  predominant  process  at 
high  pump  power  is  that  for  the  states  with  all  fields 
polarized  in  one  direction. 


3.3e  Numerical  Estimates  of  Threshold  Powers 

We  have  just  calculated  what  we  call  the  high  pump 

power  gain  which  is  the  gain  for  pump  power  well  above 

the  threshold  pump  power.  In  order  to  know  the  range 

of  pump  power  over  which  that  calculated  gain  applies 

we  must  estimate  the  threshold  power.  This  is  done  be- 

2 

low  by  finding  the  value  of  G for  which  G = Ak  . 
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— 


The  value  used  for  Ak  is  kg  since,  near  the 

II 

resonance,  kg  is  the  dominant  term  in  the  expression 

II 

for  Ak.  The  value  of  kg  is  calculated  from 


kg  = 1.71  x 10"7  z"(Cnp)  f Ag  N (3.50) 


where  f is  the  oscillator  strength  of  the  6S^  - nP 
transition.  Using  this  and  the  definition  of  G we  find 
an  expression  for  the  threshold  power. 


\-7 


1.71  x 10  ' z Us)  f As 


{2.74  x IQ-26  ~ Q2  FU_,  Cj}*5 


(3.51) 


nP 


or 


f A, 


= 1.1  x 10' 


Q2 

L S 


(3.52) 


The  values  of  ff  calculated  in  this  fashion  are  listed 
in  table  3.3  in  megawatt  units. 

Comparing  these  threshold  powers  with  the  maximum 
pump  power  of  ^ 1 Mw  we  see  that  can  be  considerably 
larger  than  threshold  power  for  10P^,  9P and  SP^.  It 
is  about  equal  to  the  threshold  for  7P^  but  is  well  be- 


low it  for  6P, 


Since  the  J = 3/2  lines  are  all  pre- 


dicted to  have  lower  gains  than  their  J = 1/2  counter- 
parts and  also  higher  thresholds,  we  consider  only  the 
J = 1/2  P states  in  the  following  discussions. 
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The  gains  calculated  above  can  be  used  for  I >>E 

P T 

but  what  is  expected  when  ^ ffT?  Since  the  gain  is 
linear  in  power  at  high  power  yet  has  a threshold,  the 
slope  of  £n  Eg  vs.  ncl’®P  must  be  9reater  near  threshold 
than  at  high  pump  powers  . 

3.4  Summary 

The  results  of  these  calculations  of  parametric 

gain  can  be  summarized  as  follows:  the  gain  profile 

versus  signal  wavelength  is  shown  in  fig.  3.9  for  the 

case  of  exact  resonance  pumping  and  has  a characteristic 

Doppler  width  centered  on  the  atomic  transitions  6S,  - nP . 

H 

The  profile  for  near-resonant  pumping  is  shown  in  fig. 
3.10.  The  gain  has  been  calculated  for  the  lOPjj  line 
for  the  experimental  conditions  encountered  and  the  gains 
of  the  other  lines  have  been  calculated  relative  to  the 
10Pi^  gain.  These  are  listed  in  table  3.3.  The  threshold 
powers,  which  have  been  estimated  assuming  negligible 
populations  in  states  other  than  6S,  are  also  listed  in 
table  3.3. 

The  following  two  chapters  will  expand  upon  the 
theory  just  presented  and  apply  it  to  the  analysis  of 
experimental  results.  Specifical1  , chapter  four  will 
deal  with  THG  and  TP A and  chapter  five  will  discuss  the 
parametric  process . Factors  such  as  saturation  of  the 
various  transitions,  Stark  shifts  induced  by  the  laser 
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radiation,  and  photoionization  from  the  9D^  level, 
which  might  modify  the  results  presented  above,  will 
be  discussed. 


CHAPTER  4 


THIRD  HARMONIC  AND  TWO 
PHOTON  ABSORPTION  EXPERIMENTS 


The  topic  of  this  chapter  is  the  description  of  a 
detailed  model  of  third  harmonic  generation  (THG)  and 
two  photon  absorption  (TPA)  and  a comparison  of  the  pre- 
dictions of  the  model  with  our  experimental  results.12 

A glance  at  some  typical  data  such  as  that  in  fig. 
4.1  reveals  that  neither  THG  nor  TPA  exhibits  the  power 
dependence  predicted  for  it  in  chapter  3.  This  failure 
may  be  understood  when  the  effect  of  the  radiation 
field  on  the  atomic  state  populations  is  considered. 

If  the  fraction  (R)  of  atoms  excited  to  the  two  photon 
resonance  state  (9D,)  is  estimated  using  R = a__._  • F • T 
where  °TPA  is  the  cross  section  calculated  in  chapter  3, 

F is  the  photon  flux  and  T is  the  pulse  duration,  we 
find  that  R may  be  as  large  as  10  for  a typical  laser 
pulse.  This  large  population  redistribution  is  in 
contradiction  to  the  assumption  of  negligible  population 
shift  implicit  in  the  perturbation  calculation  of  chap- 
ter 3.  Thus  a more  relevant  model  must  include  the 
intensity  dependent  redistribution  of  the  atomic  state 
populations . 

Such  a calculation  is  presented  and  discussed  in 
the  following  section.  Subsequent  sections  will 


Log  ( Watts ) 


Fig . 4.1 

Observed  power  laws  for  resonant  THG  (top)  and  TPA 
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compare  the  model  with  experimental  measurements  and 
discuss  the  range  of  validity  of  the  model  and  present 
our  conclusions. 

4.1  Steady  State  Saturation  Model 

The  calculations  presented  here  are  based  on  a 
steady  state  model  in  which  certain  transient  effects 
associated  with  a laser  pulse  of  finite  duration  are 
ignored.  Expressions  for  THG  and  TPA  as  functions  of 
laser  power  and  frequency  will  be  derived  and  discussed. 
Initially  only  the  population  redistribution  will  be 
considered  but  later  other  effects  such  as  Stark  shifts 
and  photoionization  will  be  included.  Throughout  the 
discussion,  the  quantities  of  primary  interest  will 
be  the  lineshapes  of  the  two  photon  resonance  and  the 
power  dependences  of  THG  and  TPA. 

The  density  matrix  formulation  of  quantum  mechanics 
is  well  suited  to  the  treatment  of  process  involving 
large  population  shifts,  and,  if  a steady  state  solu- 
tion is  assumed,  produces  simple  expressions  for  TPA 
and  THG.  Such  a calculation  is  presented  in  Appendix  A 
and  the  results  are  discussed  here.  The  steady  state 
assumption  means  that  the  population  in  each  atomic 
state  is  taken  to  be  time  independent  during  the  laser 
pulse,  because  the  pulse  duration  is  much  greater  than 
the  relaxation  times.  For  simplicity,  we  assume  that 


»SK» 


only  the  ground  state  and  the  two-photon-resonant  state 
attain  significant  populations.  In  our  case  these  are 
the  6S^  and  9D^  states  of  cesium.  For  the  moment,  loss 
of  atoms  from  the  upper  state  other  than  by  stimulated 
or  spontaneous  decay  into  the  lower  state  are  neglected. 
The  incident  radiation  is  assumed  to  be  tuneable  and  mono- 
chromatic . 

The  expressions  derived  for  the  THG  and  TPA  sig- 
nals are  (see  Appendix  A) 


TPA  « 


a + S 2T2  Aw9 


Z"(x,y) 


(4.1) 


THG  a z,2<x'y>  + T^~s  z"“U,y) 


1 _ 2 


(4.2) 


where 


2a>  - /I  + S 

9 

x = ; y = 

Aw.  T Aw 

9 2 9 


; Aoj9  = 2kT/Mc  (4.3) 


- i v,  W 


<9D^ I z |mP><mP | z | 6 S ^ > 
* (%p  - (0) 


( | | 2) 2 


(4.4) 


The  quantity  S — called  the  saturation  parameter — is 
a number  which  indicates  the  population  difference  between 
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the  two  states  9D3  and  6^  for  those  atoms  exactly 

"2  '2 

resonant  with  the  radiation.  Specifically,  for  those 
atoms 


'11 


22 


(1  + S) 


-1 


(4.5) 


where  p is  the  fraction  of  atoms  in  the  ground  state 


and 


22 


is  the  fraction  in  the  excited  state.  The 


quantities  T and  appearing  in  eqn.  4.4  are  respec- 
tively the  lifetimes  of  the  diagonal  and  off  diagonal 
terms  of  the  density  matrix  or,  in  other  words,  is  the 

lifetime  of  the  excited  state  and  T is  the  inverse  of  the 

2 


homogeneous  width  of  9D, 


It  may  easily  be 


shown  that  an  alternative  expression  for  S is 


S = 2 ‘ * F ‘ T1  (4‘S 6) 

where  a is  the  maximum  two  photon  absorption  cross 
M 

section  of  an  atom  in  the  ground  state,  is  the  life- 
time of  the  9D^  state  and  F the  photon  flux  in  number 
per  cm^  per  second. 

We  now  discuss  in  detail  the  expressions  for  THG  and 
TP A.  For  this  purpose  it  is  convenient  to  consider 
three  ranges  of  laser  power.  The  first  is  the  low  power 
range  where  S<<1  and  y<<l. 
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LOW  LASER  POWER  S<<1  y<<! 

Whan  S<<1 , the  quantity  (p^  _ p22^  near^y  unity 
according  to  eqn.  4.5,  so  there  are  few  atoms  in  the 
excited  state.  The  value  of  y is  also  much  less  than  1 if 
it  is  assumed  that  (T2  Ato^)  ^ is  small.  This  is  equiva- 
lent to  assuming  that  the  homogeneous  line  width  is  much 
less  than  the  Doppler  width.  Such  is  not  always  the  case, 
of  course,  but  under  the  conditions  of  interest  to  us  it 
is  likely  to  be.  In  either  case,  the  equations  for  TPA 
and  THG  reduce  to 


THG  * I3  -j-i—  Z2  (x,y) 
to  Ato„ 

9 


TPA  ec 


2T2  Ato9 


Z"(x,y) 


(4.7) 


(4.8) 


where 


2to  - Q ' 

9 

x = ; y = 

Acon 


T2  Aco9 


These  are  the  same  results  as  yielded  by  the  perturbation 
calculation  of  chapter  3.  This  is  expected  since  in  both 
cases  it  was  assumed  that  (p  “ P22)  = T^e  sa^ent 

features  in  this  power  range  are  quadratic  and  cubic  de- 
pendences for  TPA  and  THG  respectively  with  lineshapes 
described  by  Z"(x,y)  for  TPA  and  z2(x,y)  for  THG  where  x 


f • 
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and  y are  independent  of  I . These  lineshapes  appear 
in  figs.  4.2  and  4.3. 

INTERMEDIATE  LASER  POWER  S>>1  y<<! 

If  (T2  • Awg ) ^ is  sufficiently  small,  there  is  a 
range  of  laser  power  over  which  S is  large  yet  y is  still 
much  less  than  unity.  In  such  a case  the  equations  4.1 
and  4 . 2 become 


THG  « I3  { Z '2 (x,y)  + i Z"2 (x,y) } 


TPA  <*  /S  Z"(x,y)  a Z"(x,y) 


(4.9) 


(4.10) 


where  y = 


/S 


T2  A“9 


Recalling  that  S is  proportional  to  we  see  that 

U) 

the  TPA  lineshape  is  changed  due  to  the  larger  value  of 
y.  The  lineshape  is  still  dominated  by  Doppler  broaden- 
ing, however,  so  the  change  is  slight. 

The  behavior  of  THG  is  more  complex.  At  linecenter 
(x  = 0)  where  Z'(x,y)  = 0,  the  third  harmonic  signal  is 
linear  in  ff  while  off  resonance  where  Z ' (x ,y) >> Z" (x ,y) 

U) 

it  is  proportional  to  _ Thus  the  linecenter  is  de- 
pressed relative  to  the  wings  and  a dip  appears  there. 
This  is  illustrated  in  fig.  4.3. 

The  physical  basis  for  this  is  that  in  this  power 
range  a hole  is  burned  in  the  longitudinal  velocity 
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distribution  of  atoms — those  atoms  which  have  the  neces- 
sary velocity  along  the  direction  of  propagation  of  the 
radiation  to  be  Doppler  shifted  into  near-resonance 
with  it  are  more  strongly  saturated  than  nonresonant 
atoms.  Thus,  if  the  quantity  (p^  - p22>  is  plotted  as  a 
function  of  the  lab  frame  resonant  frequencies,  a dip 
occurs  at  the  radiation  frequency.  This  hole  has  a width 
y when  measured  in  units  of  Aujg . Since  y is  proportional 
to  /S'  or  to  , the  width  of  the  hole  is  proportional  to 
laser  power,  but  only  a small  fraction  of  all  atoms  will 
be  saturated  since  the  hole  width  is  much  less  than  the 
Doppler  width  (y<<l) . Because  the  TPA  signal  is  pro- 
portional to  p22'  will  be  proportional  to  the  hole 
width  multiplied  by  the  density  of  atoms  resonant  at 
the  radiation  frequency.  Thus,  the  TPA  signal  should  be 
proportional  to  and  the  lineshape  should  be  nearly 
the  same  as  at  lower  radiation  intensities. 

It  is  shown  in  Appendix  A that  the  contribution  to 

the  THG  signal  from  atoms  of  a given  resonant  frequency 

2 

is  proportional  to  the  value  (p-^  “ P22^  f°r  t^lose  atoms. 
The  Z"  term  is  the  contribution  from  these  saturated 
atoms  so  THG  is  expected  to  be  proportional  to  the 
square  of  the  number  of  saturated  atoms  multiplied  by 


*U  (pll  p22} 


Z"  (x,y) 


Since  (p  - p ) is  pro- 

x 1 22 


portional  to  1/(1  + S) , this  produces  the  I power  law  at 


co 


linecenter. 
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The  Z'  term  of  eqn.  4.9,  on  the  other  hand,  is  the 
contribution  from  non-resonant  atoms  and  retains  its 

0) 

power  law  so  long  as  y<<l,  which  is  to  say  as  long  as 
only  a small  percentage  of  all  atoms  is  saturated. 

HIGH  POWER  S>>1  y>>! 

As  the  radiation  intensity  increases  and  y in  turn 
increases,  eventually  the  point  is  reached  where  y>>l. 

This  means  that  virtually  all  the  atoms  will  be  saturated 
if  x<y  (x  is  the  detuning  of  the  radiation  from  linecenter 
measured  in  the  same  units  as  y) , in  which  case  the  TPA 
signal  will  be  independent  of  intensity.  The  signal  falls 
off  only  if  x becomes  greater  than  y,  so  the  TPA  linewidth 
is  y.  The  lineshape  for  a particular  value  of  y is  shown 
in  fig.  4.2.  It  is  nearly  Lorentzian  in  the  limit  of 
large  y since  the  Doppler  width  is  then  much  smaller  than 
the  homogeneous  width. 

For  THG , the  signal  near  linecenter  (x<<y)  is  pro- 
portional to  P-'*'  . This  occurs  because  all  atoms  are 

ui 

2 

saturated  and  the  signal  is  proportional  to  (p^  ~ p^  ^ * 

or  to  I--*-.  For  x>>y  the  power  law  is  still  cubic.  The 

0)  01 

resulting  high  power  THG  lineshape  is  shown  in  fig.  4.3. 

This  discussion  of  the  power  dependence  of  line- 
shapes  and  power  laws  for  TPA  and  THG  is  summarized  in 
the  set  of  figures  4.2  - 4.5.  The  lineshapes  illustrated 
in  figures  4.2  and  4.3  have  been  referred  to  above.  They 
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show  the  lineshapes  for  TPA  and  THG  for  each  of  the  three 
power  ranges. 

The  corresponding  power  laws  are  shown  in  figs.  4.4 

and  4.5.  Fig.  4.4  shows  the  dependence  of  the  TPA  signal 

on  E^  for  x = 0 and  for  x = 1.4.  At  the  linecenter 

2 1 

(x  = 0)  the  three  power  ranges  produce  the  E^ , E^ , and 

E^  slopes  as  illustrated  in  the  figure.  When  x = 1.4,  the 

signal  saturates  to  the  same  value  as  when  x = 0 and  is 

2 

proportional  to  E^  at  low  power.  The  THG  power  depen- 
dence is  shown  in  fig.  4.5.  The  x = 0 curve  has  the 
E,3  , E,1 , and  E,  ^ slopes  discussed  above.  Notice  that  de- 
tuning  from  linecenter  improves  THG  efficiency  at  high 
powers.  This  is  expected  considering  the  dip  at  line- 
center  (see  fig.  4.3).  In  fact,  maximum  efficiency  at 
high  power  is  achieved  for  x = y.  Thus,  as  E^  increases 
and  y increases,  the  driving  radiation  must  be  detuned 
in  proportion  to  E^  to  obtain  maximum  THG.  If  this  is 
done,  the  signal  will  be  proportional  to  E^  rather  than 
E^3 . Thus,  if  only  the  contribution  from  the  9D|  two- 
photon-resonant  state  is  considered,  there  is  a limiting 
efficiency  of  THG  and  increasing  the  laser  power  further 
would  not  improve  it.  Furthermore,  this  maximum  effi- 
ciency is  obtained  by  detuning  from  linecenter  by  an 
amount  proportional  to  the  laser  intensity. 

In  practice,  it  may  not  be  necessary  to  detune  the 
radiation  source  from  linecenter  to  obtain  a E^  power 
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dependence  since  the  A.  C.  Stark  effect  will  shift  the 
frequency  of  the  atomic  resonance  by  an  amount  propor- 
tional to  the  instantaneous  intensity.13  The  shift  need 
not  be  exactly  equal  to  y to  produce  a linear  intensity 
dependence.  Any  shift  which  is  proportional  to  5^  will 
produce  the  same  linear  dependence  at  sufficiently  high 
power,  although,  of  course,  it  may  produce  less  than  the 
maximum  efficiency.  The  effect  of  Stark  shifts  (for  x = 0) 
are  shown  in  fig.  4.6.  As  the  Stark  effect  shifts  the  dip 
at  linecenter  away  from  the  radiation  frequency,  the  THG 
efficiency  is  improved.  While  this  is  occuring,  the 
exponent  of  the  power  law  may  even  be  greater  than  three 
as  in  the  case  illustrated.  At  high  powers  (y>>l,  S>>1) 
the  THG  signal  becomes  linear  in  pump  intensity. 

Thus  far  in  the  analysis  of  THG  as  a function  of 
laser  intensity,  the  role  of  the  wave  vector  mismatch 
Ak  has  been  ignored.  We  have  dealt  only  with  the  in- 
fluence of  x<  assuming  Ak  to  be  constant.  However,  chang- 
ing the  population  distribution  also  affects  Ak  because 
the  linear  polarizability  of  the  sample  at  frequencies 
w and  3u>  changes  with  population.  In  addition,  the  in- 
tensity dependent  refractive  index  may  be  significant  in 
some  situations.  So  long  as  the  sample  is  only  a fraction 
of  a coherence  length,  our  discussion  remains  valid  since 
index  matching  is  not  important.  If,  however,  a sample 


of  several  coherence  lengths  is  employed  or  if  index 
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matching  by  a second  species  of  atom  or  molecule  with  Ak  of 
opposite  sign  is  attempted,  then  the  intensity  dependence 
of  Ak  can  have  a significant  influence  in  THG.  In  our 
experiments  we  have  about  one  coherence  length  of  vapor 
and  since  Ak  increases  or  decreases  by  less  than  one 
coherence  length,  the  contribution  of  Ak  is  not  as  sig- 
nificant as  the  changes  in  x> 

Another  process  which  may  be  important  in  TPA  and 
THG  is  the  loss  of  atoms  from  the  two  state  system. 

While  this  is  at  odds  with  the  restrictions  imposed  by 
the  steady  state  assumption,  it  will  certainly  be  sig- 
nificant in  some  experimental  situations.  The  loss  can 
occur  by  spontaneous  or  stimulated  decay  of  the  upper 
state  to  states  other  than  the  lower  state  or  by  colli- 
sional  energy  transfer  or  photoionization  from  either 
resonant  state.  Consider,  for  example,  single  laser 
photon  ionization  of  the  9D,  state  of  cesium.  At  suf- 
ficient radiation  intensities  the  lifetime  of  this  state 
to  photoionization  will  become  less  than  the  low  power 
lifetime.  A significant  fraction  of  atoms  which  are  pumped 
into  9D,  will  then  be  ionized.  This  will  tend  to  reduce 
the  saturation  or  increase  the  value  of  (P-^  ~ p 22^  an<^ 
consequently  improve  THG  so  long  as  the  population  frac- 
tion is  not  depleted  too  severely. 

In  addition,  the  change  in  upper  state  lifetime 
produced  by  photoionization  may  alter  the  value  of  y. 


This  lifetime  broadening  of  the  resonance  may  change 
the  optimum  detuning  and  perhaps  the  power  law  at  high 
intensities . 1 4 


The  steady  state  model  cannot  accommodate  a full 
description  of  the  system  including  loss  of  atoms  be- 
cause the  loss  is  incompatible  with  the  steady  state  as- 
sumption. Other  methods  which  can  take  these  factors 
into  account  will  be  discussed  in  a later  section. 

4.2  Comparison  of  Model  and  Experimental  Results 

In  order  to  compute  the  lineshapes  and  power  laws 
for  the  cesium  - ruby  laser  system,  it  is  necessary  to 
know  the  values  of  line  widths,  hyperfine  structure. 

Stark  shifts,  etc.  The  estimates  of  these  will  be  dis- 
cussed first  and  then  the  predictions  of  the  steady  state 
model  based  upon  these  will  be  compared  with  observations. 

A list  of  values  for  various  parameters  (at  one 
coherence  length)  is  given  below 
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Lifetime 

25 

ns 

Homogeneous  Linewidth 

> 

6.5 

Mhz 

Doppler  Linewidth 

= 

1200 

Mhz 

2 x 

Laser  Mode  Separation 

= 

250 

Mhz 

2 x 

Laser  Mode  Width 

> 

5 

Mhz 

6SJ5 

Hyperfine  Splittings 

= 

9200 

Mhz 

9D| 

Hyperfine  Splittings 

B 

10 

Mhz 

82 


825 


Stark  Shift  (65^  - 9D^) 

Photoionization  Lifetime  9D3 

% 


Mhz  x 2 (MW) 
4 ns/I (MW) 


These  are  discussed  individually  below. 

The  9D3  lifetime  was  measured  in  two  ways.  The 

first  consisted  of  estimating  the  decay  time  of  9D  from 

\ 

oscilloscope  traces  of  the  fluorescence  at  the  9D3  - 6P3 

"1  '2 

transition  frequency  (See  fig.  4.7).  Since  the  decay 
time  is  comparable  to  or  longer  than  the  laser  pulse  at 
low  cesium  densities,  the  time  constant  can  be  estimated. 
It  is  about  30  ns  at  one  coherence  length  and  about  100  ns 
at  1/3  CL  so  it  is  density  dependent  with  a lifetime  of 
'V'  30  ns  at  the  density  of  interest.  This  is  to  be  com- 
pared to  a calculated  radiative  lifetime  of  210  ns. 4 
These  estimates  were  confirmed  by  measurements  of  the 
total  fluorescence  signal  from  the  same  transition  inte- 
grated over  the  time  from  the  beginning  of  the  laser 
pulse  to  100  ns  after  the  end  of  the  pulse  as  a function 
of  cesium  density.  This  data  (shown  in  fig.  4.8)  was  fit 
by  a simple  model  including  photoionization  of  9D,, 
spontaneous  decay,  and  collision  induced  decay  at  a rate 
proportional  to  cesium  density.  The  proportionality 
constant  for  the  latter  was  chosen  to  provide  a fit  to 
the  data  and  from  it  the  collision  induced  lifetime  was 
estimated.  Assuming  a square  laser  pulse,  the  change 
in  population  of  9D|  may  be  approximated  by 


.3  CL 


< .3  CL 


Fig.  4.7 

Oscilloscope  traces  of  fluorescence  signal  at 
various  cesium  densities  (50  ns/major  div.) 
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where  p is  the  cesium  density  in  coherence  lengths , C 1 

is  the  natural  or  radiative  lifetime  of  9Di,  B- ^ the 

2 

photoionization  lifetime,  and  D the  constant  we  wish  to 
find.  When  this  equation  is  integrated  over  time  to 
find  /q  n(t)dt  which  is  proportional  to  the  TPA  signal, 
and  the  best  estimates  for  A,  B and  C are  used,  (B~‘*-  ~ 

10  ns,  C--*-  = 150  ns,  T = 20  ns)  , a value  for  D~ ^ of  25  ns 
produces  the  curve  labeled  "Fit"  in  fig.  4.8.  This  is 
insensitive  to  small  changes  in  T and  B.  The  conclusion 
is  that  the  collision  induced  lifetime  is  25  ns/pCL. 

The  value  listed  for  the  homogeneous  linewidth  is  a 
minimum  value  deduced  from  the  9Di  lifetime.  In  all 

'i 

likelihood  the  actual  value  is  considerably  larger  than 
this  due  to  line  broadening  from  Cs  - Cs  collisions. 

For  the  present  the  homogeneous  width  will  be  written 
as  3 x 6.5  Mhz  where  8 is  a number  greater  than  unity. 

The  next  width  on  the  list,  the  Doppler  width,  is 
calculated  using  a temperature  of  210°C  corresponding  to 
the  cesium  density  at  one  coherence  length. 

The  laser  mode  separation  is  the  separation  of 
adjacent  longitudinal  modes  calculated  using  ,\ v ■= 
where  L is  the  optical  length  of  the  laser  cavity, 
value  is  confirmed  by  observing  the  beat  frequer  : . 
tween  cavity  modes. 
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The  laser  mode  width  listed  is  the  transform  limited 
frequency  spread  associated  with  a 30  ns  pulse.  This  is 
a minimum  value  but  the  actual  value  should  differ  little 
from  it. 

The  hyperfine  splittings  of  66^  and  9D|  are  those 
reported  by  Svanberg,  et.  al.15  The  structure  is  dia- 
grammed in  fig.  4.9. 

The  Stark  shift  and  photoionization  rate  are  cal- 
culated in  appendices  B and  C.  The  direction  of  the 
Stark  shift  is  toward  increasing  separation  of  the  energy 


levels  6S  and  9D,  with  increasing  laser  intensity.  The 


photoionization  rate  is  for  single  photon  ionization 


of  9D|  only,  as  this  is  much  larger  than  the  ionization 


rate  of  6S^  or  multi-photon  ionization  of  9D|  for  the 


conditions  of  interest  here. 

We  now  have  enough  information  to  calculate  the 
laser  powers  at  which  S = 1 and  at  which  y = 1.  These 


powers  will  be  labeled  and  respectively.  Earlier 


we  found  that  R,  the  fraction  of  atoms  excited  by  a 


7 2 

laser  pulse,  was  unity  for  a laser  intensity  of  10  w/cm 


where  R = a • F • T.  The  value  of  T was  30  ns  and 
TPA 


a was  the  Doppler  averaged  two-photon-absorption  cross 

X r 


section.  Using  the  expression  for  the  saturation 
parameter  given  in  eqn.  4.6,S  may  be  expressed  in  terms 
of  R as  follows. 


fc*. 


Jk 


where  is  25  ns  for  1 CL  and  T is  30  ns.  The  ratio  of 
cross  sections  is  simply  the  ratio  of  the  Doppler  width 
to  the  low  power  homogeneous  width.  This  ratio  is 

appriximately  1200/(6.5  x 6)  or  185/8.  Thus  S = 

3 7 2 

(3  x 10  )/8  for  an  intensity  of  10  W/cm  . Considering 

2 

that  the  area  of  the  beam  xs  .03  cm  , this  yields  a value 
for  Is  of 


= 104  /6  watts 


(4.13) 


This  calculation  of  Ig  was  based  on  a single  mode 
laser.  Since  there  are  typically  4-6  laser  modes, 
the  power  must  be  divided  among  them.  Assuming  the 
power  is  predominantly  shared  by  only  two  or  three  modes , 
as  is  usually  the  case,  the  value  of  1 is  changed  very 

O 

little.  If  there  are  two  modes  of  equal  power,  for  -t 

example,  1 is  unchanged  since  the  two  modes  saturate 
the  atoms  resonant  at  the  frequency  midway  between  them 
just  as  effectively  as  a single  mode  with  the  same  total 

• | 

power  at  the  intermediate  frequency.  If  there  are  three 
equal  power  modes,  lg  is  increased  less  than  50%. 

.1 


As  laser  power  is  increased  above  Eg,  the  hole  width 

increases  in  proportion  to  E^  until  it  is  equal  to  the 

Doppler  width  at  E,,  = I . Thus 

a}  y 


E 


y 


Doppler 

Homogeneous 


E 


y 


1200 

6.5  x i* 


(4.13) 


= 2 x io  //$  watts 


Again  this  is  calculated  for  a single  mode  laser. 

These  calculations  of  Eg,  E^,  and  the  lineshapes  are 
now  compared  with  data  obtained  when  the  laser  was  opera- 
ting on  a single  longitudinal  mode.  The  data  shown  in 
fig.  4.10  was  obtained  by  measuring  the  TPA  and  THG 
signals  as  the  laser  was  tuned  across  the  resonance  be- 
tween 9D^  and  the  higher  energy  63^  hyperfine  component. 
The  top  pair  (Fig.  10a  and  lOd)  and  bottom  pair  (Fig. 

10c  and  10 f)  were  taken  on  alternate  laser  shots  to 
eliminate  any  drift  in  calibration  of  laser  frequency 
between  the  two.  Pair  10a  and  lOd  was  taken  at  a laser 
intensity  well  below  E ; the  second  pair  (10b  and  lOe)  at 

D 

a laser  power  of  about  105  watts  with  a collimated  beam; 
and  the  last  pair  (10c  and  lOf)  at  10^  watts  with  the 
beam  focused  at  the  center  of  the  cell  by  an  18  cm  focal 
length  lens. 
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First,  we  compare  the  data  of  fig.  10a  with  the 
lineshape  for  TPA  computed  from  the  steady  state  model. 
Various  values  of  3 are  tried  and  the  best  fit  is  ob- 
tained for  a value  of  about  50  (see  fig.  4.11a).  That 
this  is  the  best  fit  may  be  seen  most  clearly  in  the  wings 
of  the  line.  This  value  of  3 implies  a dephasing  time 
of  1 ns  for  the  state  9D^  which  is  not  unreasonable  at 
this  cesium  density. 

Now  using  this  value  for  3,  we  fit  a curve  to  the 
data  in  fig.  10b.  The  value  of  the  saturation  parameter 
is  varied  to  obtain  the  best  fit  (shown  in  fig.  4.11b)  for 
S = 40.  This  implies  that  the  laser  power  used  to  gener- 
ate fig.  10b  should  have  been  /S  • lg  or  4 . 5 * 10^  watts 

which  compares  favorably  with  the  measured  laser  power  of 
5 

10  watts,  so  thus  far  the  steady  state  model  has  pro- 
vided satisfactory  explanations  of  the  TPA  data  at  low 
and  intermediate  powers. 

Next,  the  data  of  fig.  10c  is  considered.  The 
lineshape  is  fit  fairly  well  by  the  same  curve  as  was 
fig.  10b,  but  there  is  a shift  to  the  left  relative  to 
10a  by  300  Mhz.  This  is  in  the  direction  expected  for 
the  Stark  shift  and  the  size  is  close  to  the  375  Mhz 
calculated  for  4.5  * 105  watts  laser  power  inferred  from 
the  lineshape.  Note  that  a shift  in  10b  should  be  of  the 
same  size  since  the  same  intensity  is  used  to  fit  the 
lineshape,  but  possible  drifts  in  laser  frequency 


is  computed  without  Stark  shift;  that  in  d is  with. 


k . 4 4 


3 

3 
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Why  does  the  TP A signal  for  a focussed  beam  have  the 
same  resonance  profile  as  the  collimated  beam  of  the  same 
power  when  the  intensity  near  the  focus  (where  the  fluo- 
rescence is  observed)  must  be  much  greater  than  for  the 
collimated  beam  and  consequently  the  line  should  be 

broadened?  The  explanation  is  that  at  intensities  of 
5 2 

5 k 10  watts/cm  the  photoionization  lifetime  becomes 
shorter  than  the  laser  pulse  width  and  most  cesium  atoms 
are  ionized.  This  depletes  the  number  of  fluorescing 
atoms  and  the  TPA  signal  is  expected  to  decrease  with  in- 
creasing power.  When  the  laser  power  becomes  great  enough 
for  this  to  happen  at  the  beam  center,  the  signal  from  the 
weaker  parts  of  the  beam  dominate  and  instead  of  decreas- 
ing with  increasing  power,  the  TPA  signal  is  independent 
of  laser  power.  This  can  be  shown  by  integrating  the 
signal  across  the  beam  profile.  As  the  power  increases 
further,  the  signal  always  comes  predominantly  from  that 

part  of  the  beam  where  the  intensity  is  about  3 * 10^ 

2 

watts/cm  which  corresponds  to  the  maximum  intensity  at 
10^  watts.  Thus,  the  lineshape  should  not  continue  to 

0) 

broaden  nor  should  the  Stark  shift  increase  with  intensity 
if  a Gaussian  beam  is  used.  In  a sense  this  verifies  that 
the  calculated  value  of  the  photoionization  rate  is  of  the 
right  magnitude. 


l 
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The  conclusion  to  be  drawn  from  the  comparison  of 
TPA  lineshapes  at  various  powers  with  the  steady  state 
model,  is  that  the  agreement  is  well  within  the  bounds  of 
uncertainty  of  calculations  and  measurements  in  all  de- 
tails if  0 is  50. 

Next  the  THG  data  in  fig.  4.10e  and  4.10f  are  com- 
pared with  computed  lineshapes  for  various  values  of  S. 

It  is  found  that  S must  be  less  than  5 to  provide  a rea- 
sonable fit  to  the  data  of  lOe  (see  fig.  4.12a).  The 
value  of  40  which  was  found  to  provide  the  best  fit  for 
TPA  at  the  same  power  produced  the  result  in  fig.  4.12b. 
Thus  there  is  a factor  of  10  or  so  difference  in  the 
saturation  parameter  deduced  from  the  TPA  and  THG  line- 
shapes  at  this  power.  Another  notable  difference  is  that 
there  is  substantially  more  scatter  in  the  THG  data  than 
in  TPA  data. 

The  scatter  in  fig.  4.l0f  is  likewise  large  and 
makes  it  almost  impossible  to  determine  the.  value  of  S 
which  provides  the  best  fit  to  the  data.  It  is  clear, 
however,  that  the  resonance  is  considerably  broadened  or 


shifted  for  this  focussed  beam  case. 


With  the  laser  operating  on  a single  longitudinal 
mode,  the  available  power  was  insufficient  for  measure- 
ments of  the  dependence  of  THG  on  laser  power.  The  higher 
power  available  from  multimode  operation  makes  this  pos- 
sible and  also  provides  information  about  lineshapes  at 


Fig.  4.12 

THG  data  fit  by  steady  state  model.  Lineshape  in  a 
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higher  intensities.  This  data  is  reproduced  in  figs. 

4.13  and  4.14.  The  laser  ; owers  in  fig.  4.13  are  approxi- 

4 ^ 

mately  10  watts  for  the  upper  pair  (a,e) , 10  watts  for 

the  second  pair,  and  106  watts  for  the  third  pair  all  with 

a collimated  beam.  The  final  pair  (d,h)  is  for  10^  watts 

and  a focussed  beam.  Fig.  4.14a  presents  data  on  THG  as 

a function  of  laser  power  when  the  laser  is  tuned  about 

1 cm-"*'  from  the  two  photon  resonance,  14. b is  the  same 

measurement  except  the  laser  is  tuned  to  the  linecenter 

of  the  resonance  between  the  upper  ground  state  hyperfine 

level  and  9D and  14. c is  the  TPA  signal  as  a function 

of  laser  power  at  the  same  laser  frequency  as  14. b. 

The  data  of  fig.  4.13  has  been  compared  to  calculated 
lineshapes  much  as  in  the  single  mode  case  except  three 
modes  of  equal  power  separated  by  700  Mhz  have  been  used 
to  approximate  the  actual  multimode  laser  spectrum.  The 
characteristics  of  the  computed  lineshapes  are  similar 
to  those  for  the  single  mode  case.  At  10^  watts  (figs. 
4.13b  and  f)  there  is,  as  before,  noticeable  broadening 
of  the  TPA  lineshape  compared  to  the  lower  power  case, 
and  the  saturation  parameter  is  found  to  be  about  20 
(see  fig.  4.15b)  which  is  in  good  agreement  with  the 
value  of  40  estimated  for  the  same  power  in  the  single 
mode  case.  Once  again,  the  THG  lineshape  shows  less  evi- 
dence of  saturation  than  TPA  and  the  saturation  parameter 
deduced  from  the  THG  lineshape  cannot  be  much  greater  than 
unity. 
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Fig.  4.13 

TPA  (left)  and  THG  (right)  data  for  multimode 
laser  at  10^  watts  (a,e)  , 10^  watts  (b,f), 

106  watts  (c,g) , and  106  watts  focussed  (d,h) 
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At  106  watts  (figs.  4.13c  and  g) — a power  which  was 
not  possible  in  single  mode  operation — there  is  evidence 
of  further  broadening  of  TPA  and  especially  of  THG.  A 
saturation  parameter  of  140  provides  a fit  to  TPA  (see 
fig.  4.15c)  but  a smaller  value  of  perhaps  40  yields  a 
more  satisfactory  match  to  the  THG  data  (fig.  4.15g). 

Just  as  in  the  single  mode  case,  the  focussed  beam 
lineshape  (fig.  4.13d  and  h)  is  the  same  as  the  unfocussed 
shape  for  TPA  but  again  THG  exhibits  considerable  broad- 
ening. 

Next  consider  the  dependence  of  the  TPA  and  THG  sig- 
nals on  laser  power.  The  data  is  shown  in  fig.  4.14. 

In  fig.  4.14a  it  is  evident  that  THG  is  proportional  to 
the  cube  of  the  laser  power  off  resonance  while  fig. 

4.14b  indicates  a near  quadratic  power  law  on  linecenter. 
The  dependence  of  TPA  on  laser  power  (fig.  4.14c)  is  seen 

5 

to  be  linear  at  powers  below  5 * 10  watts  and  less  than 
linear  at  higher  powers. 

These  are  to  be  compared  with  the  calculated  power 
laws  for  THG  and  TPA  shown  in  figs.  4.16  and  4.17,  which 
are  computed  for  linecenter  tuning  corresponding  to  fig. 
4.14b  and  c,  and  are  based  on  the  steady  state  saturation 
model  using  the  three  laser  modes  and  including  hyperfine 
structure  and  Stark  shifts  with  & = 50  but  not  including 
photoionization  or  other  loss  of  atoms  from  9D^. 


Fig.  4.16 

TPA  power  law  computed  with  (a)  and  without  (b)  A.C.  Stark  shift 


THG  power  law  computed  with  (a)  and  without  (b)  A.C.  Stark  shift 
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The  THG  curve  of  fig.  4.17  exhibits  a cubic  power 
law  below  10^  watts  and  a quadratic  dependence  from  10~* 
to  106  watts.  This  is  in  good  agreement  with  the  ob- 
served behavior.  You  may  have  noticed  that  this  com- 
puted power  dependence  is  quite  different  from  that  pre- 
sented earlier  when  discussing  the  characteristics  of 
THG  in  the  steady  state  approximation  (see  fig.  4.5).  The 
reason  is  that  fig.  4.5  was  calculated  assuming  a single 
laser  mode,  a single  ground  state  hyperfine  level,  and  a 
value  for  6 of  unity  so  it  is  a somewhat  unrealistic 
model  used  to  dramatize  the  effects  of  saturation  and 
line  broadening. 

The  calculated  TPA  power  dependence  shown  in  fig. 

5 

4.16  is  quadratic  below  10  watts  and  approximately 

5 6 

linear  in  the  power  range  10  - 10  watts.  The  measured 

5 

law,  in  contrast,  is  linear  below  10  watts  and  approxi- 

H 5 6 

mately  proportional  to  in  the  range  10  - 10  watts 

(see  fig.  4.14c).  Additional  data  indicates  a quadratic 

3 

dependence  below  10  watts  or  so. 

We  can  summarize  the  comparison  between  our  measure- 
ments and  the  steady  state  saturation  model  as  follows. 

The  TPA  lineshapes  fit  the  model  well  for  both  the  single 
and  multimode  cases  if  3 is  50.  The  lineshapes  of  THG 
in  general  show  more  scatter  than  the  TPA  data  and  it  is 
difficult  to  match  each  curve  with  a calculated  one. 

It  is  clear,  nevertheless,  that  the  resonance  broadens 


with  increasing  power  and  that  the  saturation  parameter 
necessary  to  fit  THG  data  is  lower  than  for  TPA  at  the 
same  power.  The  dependence  of  THG  on  laser  power  also  is 
in  agreement  with  the  calculated  value,  being  quadratic 
at  linecenter  and  cubic  off.  The  TPA  power  law,  on  the 

other  hand,  is  puzzling.  It  is  expected  to  be  quadratic 
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below  10  watts  but  instead  is  linear  from  10  to  10 

watts . Once  again , it  should  be  pointed  out  that  the 

calculations  based  on  the  steady  state  saturation  model 

do  not  take  account  of  loss  of  atoms  from  the  two  state 

system  nor  do  they  consider  transient  effects  which  are 

important  at  sufficiently  short  times.  Perhaps  these 

effects  could  account  for  the  discrepancies.  The  role 

of  each  in  TPA  and  THG  is  discussed  below. 

Consider  first  the  transient  behavior  which  will  be 
important  if  the  pulse  duration  is  sufficiently  short. 

If  the  laser  pulse  is  approximated  as  a square  pulse, 
then  for  a time  after  it  is  turned  on  the  atomic  system 
will  exhibit  precession  in  which  the  population  oscillates 
between  the  ground  state  and  the  two-photon-resonant  state 
at  the  frequency 


= |^S/(2T1T2)j: 


(4.14) 


for  those  atoms  resonant  at  the  laser  frequency.  Those 
atoms  which  are  not  resonant  will  precess  with  a higher 
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frequency  but  lower  amplitude.  For  = 25  ns,  T2  = 1 ns 
(corresponding  to  3 = 50)  , and  S = 1,  the  value  of  w is 
(7  ns)  ^ , so  the  precession  frequency  is 


(4.15) 


This  precession  is  damped  by  dephasing  the  atoms 
which  occurs  on  a characteristic  time  T2  and  by  decay 
from  9Dj  which  occurs  on  a time  scale  Ti.  Only  after  the 
precession  is  damped  and  the  populations  have  reached 
their  equilibrium  values  is  the  system  described  by  the 
steady  state  model.  How  long  does  this  take  and  can 
transient  behavior  be  of  significance  in  our  experimental 
situation?  At  times  greater  than  T2  (1  ns) , after  the 
precession  is  damped  out,  simple  rate  equations  apply. 

The  time  required  to  reach  equilibrium  is  then  approxi- 
mately T^/(l  + S)  or  Tj  whichever  is  greater.  For  those 
atoms  which  are  highly  saturated  at  equilibrium  the 
precession  is  significant  for  a nanosecond  or  so  after 
which  the  equilibrium  value  is  achieved.  For  those  atoms 
which  are  saturated  little  at  equilibrium,  the  population 
difference  may  not  reach  its  equilibrium  value  until  a 
time  of  the  order  of  T^  after  the  leading  edge  of  the 
laser  pulse. 

The  TPA  signal  will  be  influenced  little  by  preces- 


Li 

II 


sion  because  it  is  so  short  in  duration  but  the  slow 


build  up  of  population  in  9D^  for  those  atoms  which  are 
only  slightly  saturated  will  tend  to  decrease  the  TPA 
from  the  value  predicted  by  the  steady  state  model.  The 
effect  is  not  of  major  significance,  however,  and  does 
not  explain  the  linear  TPA  power  dependence.  The  THG 
signal  may  be  more  significantly  affected.  It  is  possi- 
ble that  the  signal  generated  during  the  first  nanosecond 
or  so  when  precession  is  most  significant  could  be  as 
large  as  that  generated  during  the  much  longer  remainder 
of  the  pulse  since  the  signal  then  may  be  strongly  satu- 
rated. This  would  tend  to  make  the  dip  at  linecenter  less 
pronounced.  On  the  other  hand,  the  slow  build  up  to 
equilibrium  will  tend  to  increase  the  THG  coefficient 

off  resonance  since  the  value  of  (p,  - p ) will  on  the 

11  22 

average  be  greater  than  its  equilibrium  value,  so  it  is 
difficult  to  assess  the  role  of  transient  behavior  on 
THG  without  resorting  to  numerical  calculations  which  we 
do  not  attempt  here. 

Consider  next  the  loss  of  atoms  from  the  two  state 
system.  This  loss  may  occur  in  several  ways.  Photoioniza- 
tion of  9D^  by  the  laser  radiation  field  will  be  signifi- 

5 

cant  at  laser  powers  above  % 10  watts  as  was  discussed 
previously  in  connection  with  the  TPA  lineshape.  In 
addition  stimulated  and  spontaneous  decay  of  9D^  atoms  to 
any  of  several  lower  energy  P and  F states  is  possible  and 
several  parametric  processes  are  known  to  occur.  The  loss 
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of  atoms  will  affect  the  TPA  signal,  Ak,  and  in  dif- 

ferent  ways.  It  is  clear  that  TPA  will  be  decreased  by 
the  loss  of  9D^  atoms  but  the  response  of  x and  Ak  are 
more  difficult  to  ascertain.  For  example,  x will  under 
some  circumstances  be  increased  since  the  loss  of  9D| 
atoms  tends  to  increase  the  value  of  (P-^  “ p22^  w^ile  it 
will  be  decreased  if  the  loss  of  atoms  from  6S  is  severe 
enough  to  decrease  the  population  difference.  Of  course, 
the  value  of  Ak  depends  on  the  details  of  the  population 
redistribution  and  like  x is  time  dependent.  It  seems 
likely  that  the  loss  of  atoms  from  the  two  state  system 
and  the  resulting  departure  from  the  steady  state  approxi- 
mation can  account  for  the  large  scatter  in  THG  data  com- 
pared with  the  TPA  data.  The  losses  affect  THG  in  a much 
more  complex  manner  than  TPA  and  there  may  be  considerable 
differences  in  the  time  development  from  shot  to  shot 
considering  the  variety  of  loss  mechanisms. 

Of  course,  it  is  conceptually  possible  to  include  the 
loss  of  atoms  due  to  photoionization  and  other  processes 
plus  Stark  shifts  in  the  equations  of  motion  of  the 
density  matrix  and  to  solve  them  for  the  instantaneous  THG 
and  TPA  coefficients  and  for  Ak  and  then  to  integrate  over 
time  to  find  the  TPA  and  THG  signals.  This  procedure  re- 
quires extensive  application  of  numerical  techniques  in 
solving  and  integrating  the  equations  with  little  insight 
to  be  gained  about  the  roles  of  various  processes.  In 
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addition,  it  is  difficult  to  model  the  loss  of  atoms 
accurately  because  of  the  complexity  of  the  system  in 
which  decay  from  9D^  to  several  lower  energy  P and  F 
states  is  possible  and  several  parametric  processes  occur. 


4.3  Conclusions 


Because  of  the  difficulties  in  characterizing  all  of 
the  processes  which  play  a role  in  the  production  of  TPA 
and  THG  in  our  experimental  situation,  we  have  chosen  to 
identify  those  physical  processes  which  are  of  primary 
significance  and  have  constructed  a simple  model  includ- 
ing only  them.  This  model  is  quite  successful  in  account- 
ing for  our  experimental  observations.  Saturation  broad- 
ening of  the  two  photon  resonance  for  TPA  and  THG  with 
laser  power  have  been  demonstrated  and  are  in  good  agree- 
ment with  the  predicted  values  over  a range  of  33  x io^ 
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to  33  x 10  watts/cm  . The  broadening  of  TPA  is  somewhat 
greater  than  for  THG  at  any  given  power  but  this  could  be 
due  to  loss  of  atoms  from  the  two  level  system.  In 
addition,  the  THG  power  law  fits  the  model  both  on  and  off 
resonance  over  the  power  range  examined  and  the  Stark 

shift  is  verified.  Only  the  TPA  power  law  does  not  fit  the 

3 2 7 2 

model.  From  33  x 10  watts/cm  to  33  x 10  watts/cm  the 
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dependence  is  linear  in  laser  power  rather  than  quadratic 
as  expected.  This  we  are  unable  to  account  for. 


CHAPTER  5 


PARAMETRIC  FREQUENCY  CONVERSION 

The  theory  of  parametric  frequency  conversion  was 
outlined  in  chapter  three.  Briefly,  the  process  may  be 
considered  to  be  two  coupled  three-wave  mixing  processes 
which  can,  under  certain  conditions,  result  in  an  expo- 
nential growth  of  the  intensity  of  the  signal  and  idler 
waves  at  the  expense  of  the  pump  wave.  The  signal  and 
idler  typically  have  frequencies  different  from  the 
pump  and,  before  being  amplified,  much  lower  intensities 
than  the  pump. 

Several  authors  have  reported  observations  which 
they  attribute  to  parametric  frequency  conversion,  in 
particular  in  vapors  of  rubidium,  potassium,  and  water 
pumped  by  high  powered  lasers.16  In  the  case  of  cesium 
vapor  pumped  by  a ruby  laser,  Leung  observed  radiation 
at  two  signal  frequencies  emanating  from  the  sample  when 
only  the  pump  (laser)  was  incident.4  We  have  extended 
Leung's  observations  to  include  detailed  observations  of 
four  signal  and  two  idler  waves. 

The  calculations  presented  in  chapter  three  indicate 
that  under  conditions  achieved  experimentally  (1  CL  of 
Cs  and  1 MW  at  693.5  nm)  the  conversion  of  energy  from 
pump  to  signal  and  idler  should  approach  100%.  This  is 
not  achieved  in  practice.  To  understand  why  not  and  to 
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discover  whether  the  observed  signals  are  indeed  pro- 
duced by  a parametric  process,  we  have  made  a series  of 
measurements  of  signal  and  idler  intensities  as  func- 
tions of  laser  intensity  and  tuning,  and  of  cesium  den- 
sity. The  temporal  characteristics  of  the  signal  has 
also  been  studied.  To  our  knowledge,  this  is  the  first 
systematic  study  of  a parametric  process  in  a simple 
atomic  system. 

The  following  section  presents  most  of  our  parametric 
data.  In  succeeding  sections  we  will  compare  this  data 
with  the  theory  of  chapter  three  and  discuss  the  modi- 
fications of  that  theory  necessitated  by  population  re- 
distribution and  other  effects  not  considered  in  the 
simple  perturbation  calculation.  Finally,  we  examine 
the  question  of  whether  we  have  observed  a parametric 
process  of  the  type  postulated  and  conclude  that  we  have. 

5.1  The  Data 

When  the  laser,  at  a power  of  1 MW,  was  tuned  to 
the  6Sjj  - 9Dj|  two-photon  resonance  of  Cs  and  passed 
through  approximately  one  coherence  length  of  vapor, 
strong  radiation  near  the  wavelengths  348  nm,  362  nm, 

389  nm,  and  457  nm  was  observed  emanating  from  the  cell 
colinear  with  the  laser  radiation.  These  are  the  wave- 
lengths associated  with  transitions  from  nP  to  6S  when 
n = 10,  9,  8,  and  7 respectively.  We  shall  refer  to 
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these  as  the  signal  waves.  The  wavelengths  were  not 
measured  precisely  enough  to  resolve  the  fine  structure 
of  the  P states  so  the  signal  waves  are  indicated  in  fig. 
5.1  without  specifying  which  j level  is  involved. 

The  maximum  combined  power  of  the  four  signal  waves 
was  about  ten  watts.  Each  wave  was  collimated  and  co- 
linear  with  the  pump  beam  and  had  nearly  the  same  dia- 
meter as  the  pump  beam.  No  strong  signals  other  than 
those  mentioned  were  present  in  the  range  of  350  to  450 
nm. 

Idler  waves  at  1.4  pm  and  3 ym  corresponding  to 
transitions  from  9D^  to  7P  and  to  8P  were  also  studied. 
Again,  fine  structure  of  the  P state  was  unresolved.  The 
other  idlers  expected  at  the  9D^  - 9P  and  9D^  - 10P  trans- 
ition frequencies  were  outside  the  spectral  range  of  our 
detector  and  so  were  not  observed. 

The  temporal  behavior  of  the  total  combined  signal 
wave  power  is  illustrated  in  fig.  5.2.  The  first  pulse 
of  each  trace  is  the  parametric  signal  and  the  second  is 
the  pump  pulse  delayed  by  about  35  nanoseconds.  They 
were  recorded  with  a photodiode  - oscilloscope  combination 
with  a subnanosecond  risetime.  Notice  that  the  signal 
wave  pulses  follow  quite  closely  the  pump  pulses,  but 
the  modulation  of  the  pump  is  generally  exaggerated  in 
the  signal.  Further  analysis  of  similar  data  reveals  a 
power  dependent  shift  in  the  time  delay  between  the  points 
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of  maximum  power  of  the  two  pulses.  The  delay  increases 
by  about  two  to  four  nanoseconds  when  laser  power  is 
doubled  in  the  range  of  a megawatt,  but  the  pulse  shape 
is  not  substantially  altered. 

The  strength  of  the  individual  signal  waves  were 
measured  as  functions  of  laser  (pump)  power,  cesium 
density,  and  laser  frequency.  This  data  is  presented 
in  figures  5.3  - 5.6.  Figure  5.3  depicts  the  dependence 
of  the  10P  signal  (348  nm)  energy  on  laser  or  pump  power 
for  three  values  of  cesium  density.  Notice  that  the 
energy  ranges  over  more  than  nine  decades  which  is 
typical  for  the  signal  waves;  the  lower  limit  of  detect- 
able signal  being  about  ten  photons  and  the  maximum  power 
produced  being  typically  about  a few  watts.  Similar  data 
for  the  other  three  signal  waves  are  shown  in  fig.  5.4. 
Fig.  5.5  shows  the  density  dependencies  of  the  8P  and  10P 
signal  waves.  Like  the  power  function  curves,  these  indi- 
cate that  the  signal  gain  is  large  at  low  signal  powers 
but  saturates  at  a signal  power  of  a few  watts.  The  de- 
pendence on  laser  frequency  of  the  10P  signal  may  be  seen 
in  fig.  5.6.  The  upper  curve  is  for  a laser  power  of 
1 MW;  the  lower  is  for  100  kW.  Saturation  is  again  evi- 
dent in  flat  tops  of  these  curves. 


5.2  Comparison  of  Data  and  Theory 

From  the  data  presented  in  figures  5.3  - 5.5  we 
deduce  the  small  signal  gains  for  the  four  parametric 
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Fig.  5.3 

10P  parametric  signal  as  function  of  pump  power.  Note  scale  change  at  200  kW. 
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lines  and  estimate  the  pump  threshold  powers  and  these 
values  are  then  compared  with  those  computed  in  chapter 
three.  From  figs.  5.3  and  5.4,  the  small  signal  gains  at 
a density  corresponding  to  one  coherence  length  are  found 
to  be  in  proportion  of  1,  .5,  .3,  .2  and  for  the  10P , 9P, 
8P  and  7P  lines  respectively  and  the  gain  for  the  10P  line 
is  approximately  200  decades/MW.  From  fig.  5.5,  the 
experimental  small  signal  gain  as  a function  of  density  at 
a constant  laser  power  of  one  megawatt  is  750  decades/CL 
for  the  10P  line  and  250  decades/CL  for  the  8P  line.  The 
experimental  parametric  thresholds  are  estimated  from 
figures  5.3  and  5.4  to  lie  in  the  range  of  10  - 20  kW  for 
all  four  lines.  These  measured  values  are  listed  in  Table 
5.1  next  to  the  values  calculated  in  chapter  three. 

Recall  that  the  calculated  expression  for  the  gain 
of  the  10P^  line  was 

gaine  = 51  • I • N (5.1) 

where  1 is  in  Mw  and  N is  in  coherence  lengths,  so  the 
small  signal  power  gain  is  22  decades/MW  if  N is  unity 
or  22  decades/CL  if  1 is  one  megawatt.  The  gains  of  the 
other  lines  were  calculated  relative  to  this  and  since  the 
j = h lines  always  have  greater  gain  than  the  j = -I  lines, 
only  the  former  are  listed.  Threshold  powers  were  esti- 
mated by  calculating  the  power  at  which  parametric  gain 
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overcomes  linear  loss  of  the  signal  due  to  6S  - nP  ab- 
sorption. 

Comparison  of  the  power  gains  listed  in  the  table 
indicates  that  the  experimental  small  signal  power  gain 
of  10P  is  a factor  of  ten  larger  than  the  calculated 
gain  and  the  gain  measured  for  the  other  three  signals 
decrease  with  decreasing  signal  frequency  while  the  cal- 
culations indicate  the  gain  should  increase.  The  factor 
of  10  difference  in  absolute  gains  is  within  uncertain- 
ties of  calibration  or  calculation,  but  the  trends  of  gain 
with  signal  frequency  are  significant  because  they  are 
independent  of  the  uncertainties. 

Comparison  of  the  density  gains  reveals  a similar 
pattern.  The  ratio  of  the  measured  10P  to  8P  gain 
is  about  3:1  while  the  calculated  ratio  is  about  3:4.  In 
this  instance,  however,  the  difference  of  750  to  22  in 
absolute  gain  of  10P  is  even  larger  them  for  the  power 
gain.  Threshold  powers,  listed  in  the  last  two  columns, 
also  show  discrepancy  between  theory  and  experimental 
values. 

Finally,  perhaps  the  most  striking  feature  of  the 
data  in  figures  5.3  to  5.5,  is  the  saturation  in  signal 
power  evident  in  each  instance.  Saturation  always  occurs 
at  powers  of  about  a watt  although  the  shape  is  somewhat 
different  for  the  density  curves  (fig.  5.5)  them  for  the 
power  curves  (figs.  5.3  and  5.4). 


10P 
1 MW 


3 .2  .1  0 .1  .2  -3 

cm’1 


Fig.  5.6 

Variation  of  10P  signal  intensity  with 
laser  frequency  at  105  watt3  (lower)  and 
106  watts  (upper) . 


In  calculating  the  parametric  gain,  lineshape  and 


threshold  power  in  chapter  three,  we  assumed  that  all 
cesium  atoms  were  in  the  ground  state.  However,  since  we 
have  shown  that  population  in  other  states,  especially 
9D^,  is  important  in  other  nonlinear  processes,  we  should 
modify  the  theory  to  include  that  possibility.  When  9D^ 
is  populated  as  well  as  6S,  the  susceptibilities  for  sig- 
nal and  idler  frequencies  become 


Xc  = n6Xs(6S)  + NqXc(9D) 


9 AS 


(5.1) 


XI  = N6XI(6S)  + N9xi(9D) 


where  Ng  is  the  number  of  atoms  in  6S  and  Ng  is  the  number 
in  9D^.  Thus2 


N. 


N, 


® (^96  — 2u>p)  (fipg  ” (^gg  — 2u)p)  (^gp  — Wj) 


N, 


N, 


1 ^96  “ 2u)P^  ^P6  ~ ^9P  “ ^96  “ 2u)P^ 

(5.2) 

* 

The  Doppler  averaged  susceptibilities  <XS>  and  <XI>  each 
have  a narrow,  tall  spike  at  = wp{.  when  2wp  = 
similar  to  that  described  and  graphed  in  chapter  3.  This 
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spike  comes  from  the  9D  term  of  <XS>  and  the  6S  term  of 

it  it 

<X1>‘  The  product  <XS><XI>  therefore  also  has  a large 

spike  at  wg  = wpg  of  amplitude  <xs  (9D)  ><Xj  (6S)  > • Ng  • Ng 

which  causes  a large  gain  at  linecenter.  This 

replaces  the  dip  which  appears  at  linecenter 

when  Ng  or  Ng  is  zero.  The  maximum  amplitude  of  the  gain 

is  then  larger  than  the  previously  calculated  maximum  by  a 

factor  of  approximately  (Aw  , /Aw  ) • ZnTnT 

Doppler  Homogeneous  6 9 

so  the  gain  is  very  sensitive  to  the  homogeneous  widths 

and  hence  to  laser  power  and  cesium  density. 

The  complete  description  of  the  gain  profile  and  its 

variation  with  power  and  density  requires  more  detailed 

knowledge  of  the  homogeneous  widths  than  is  available.  We 

can,  however,  by  making  certain  reasonable  assumptions, 

determine  whether  including  the  9D|  population  leads  to  a 

more  accurate  description  of  the  parametric  process.  In 

the  range  of  laser  power  over  which  the  parametric  signals 
4 6 

are  observed  (10  - 10  watts) , we  found  in  analyzing  THG 

and  TPA  that  the  fraction  of  atoms  in  9Da  varies  from  per- 
haps  .1  to  .5.  Depending  on  the  linewidths  involved,  even 
the  lower  value  of  .1  could  be  sufficient  to  significantly 
enhance  the  parametric  gain.  In  this  case,  if  the  homo- 
geneous linewidths  of  the  P states  increase  with  decreas- 
ing principal  quantum  number,  the  gain  would  tend  to  be 
larger  for  higher  frequency  signal  waves  as  observed  (see 
table  5.1)  instead  of  lower  as  would  be  the  case  when 

Ng  = 0. 


The  effect  of  9D^  population  also  provides  a reason- 
able explanation  of  the  difference  in  the  power  and  densi- 
ty gains  (see  table  5.1).  The  density  gain  is  measured 
at  low  densities  and  high  power  while  the  power  gain  is 
measured  at  high  density  and  low  power.  The  density  gain 
should  *hus  be  greater  than  the  power  gain  both  because 
the  homogeneous  widths  are  less  at  low  density  and  be- 
cause the  9D  population  is  larger  at  high  laser  intensi- 
ties . 

When  9D  is  populated,  the  increase  in  parametric 
gain,  the  decrease  in  linear  absorption  of  the  signal 
wave  and  the  contribution  to  the  idler  linear  gain  all 
contribute  in  lowering  threshold  powers  from  those  calcu- 
lated for  a 6S  only  population  (see  table  5.1)  so  it  is 
not  surprising  that  the  thresholds  are  lower  than  had 
been  predicted. 

We  conclude  that  including  a 9D|  population  in  the 
simple  perturbation  calculation  of  parametric  gain  results 
in  a significantly  improved  description  of  the  observed 
low  power  gains  and  threshold  powers.  Exact  calcula- 
tions are  not  possible  given  our  lack  of  knowledge  of 
homogeneous  linewidths  but  reasonable  values  suffice. 

Might  this  or  a similar  mechanism  also  explain  the 
saturation  of  the  signals  at  about  a watt?  If  so,  judging 
from  the  shape  of  the  gain  curves  (see  figs.  5.3  - 5.6),  it 
must  depend  primarily  on  the  strength  of  the  signal  or 
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idler  itself  rather  than  the  laser  power  or  the  cesium 
density.  This  suggests  that  population  of  the  P state, 
either  by  absorption  of  the  signal  wave  by  atoms  in  6S, 
or  by  emission  of  idler  by  9D^  atoms,  is  responsible. 

Such  a P state  population  might  reduce  the  gain  substan- 
tially. Under  certain  conditions  the  gain  can  become 
zero  when  the  populations  of  6S,  9D|,  and  nP  are  equalized. 
Calculations  show  that  approximate  equalization  of  popu- 
lations can  occur  via  signal  wave  absorption  by  6S  atoms 
if  the  signal  powers  are  approximately  1,  15,  70,  180  watts 
for  the  7P,  8P,  9P,  or  10P  lines  respectively.  Although 
180  watts  is  considerably  larger  than  the  actual  saturated 
value  of  1 - 10  watts,  this  is  a more  likely  saturation 
mechanism  than  population  of  nP  by  decay  of  9D^  since 
population  of  7P  by  9D|  - 7P  transitions  would  require 
the  emission  of  500  watts  of  3 ym  radiation  whereas  we 
measured  about  .1  watts  at  this  wavelength. 

Thus  far  we  have  compared  our  observations  only  with 
calculations  of  parametric  gain.  Other  processes  such  as 
lasing  or  Raman  gain  at  the  idler  frequency  in  conjunc- 
tion with  difference  frequency  generation  at  the  signal 
frequency,  or  lasing  at  the  signal  frequency  must  be  con- 
sidered. These  are  compared  with  the  parametric  process 
below  and  rejected  as  the  dominant  process. 

The  existence  of  populations  in  excited  atomic 
states  admits  the  possibility  of  lasing  action  at  various 
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frequencies  including  idler  or  signal.  The  expressions 
for  gain  described  in  chapter  three  take  into  account 
any  changes  in  linear  loss  or  gain  due  to  population 
shifts  via  the  inclusion  of  Ak  so  lasing  should  not  be 
considered  a process  independent  of  parametric  mixing. 
Nevertheless,  if  lasing  gain  is  much  larger  than  para- 
metric gain,  the  process  resembles  lasing  rather  than 
parametric  gain.  In  this  extreme,  gain  in  the  reverse 
direction  would  be  equal  to  that  in  the  forward  direction, 
for  example,  and  phase  mismatch  would  lead  to  modulation 
of  waves  produced  by  mixing  the  lasing  and  the  pump  waves. 
Raman  gain  shares  those  characteristics.  Parametric  gain, 
in  contrast,  is  much  larger  in  the  forward  than  reverse 
direction  and  modulation  is  absent. 

To  verify  that  we  are  observing  a parametric  process 
rather  than  lasing  or  Raman  gain,  we  have  looked  for  sig- 
nal and  idler  waves  propagating  in  the  reverse  direction 
and  find  that  the  combined  signal  wave  intensity  is  less 
than  a thousandth  that  in  the  forward  direction.  Thus 
lasing  at  the  signal  frequency  is  not  significant.  The 
ratio  for  the  3 ym  idler,  on  the  other  hand,  is  between 
3 and  10  although  there  is  considerable  uncertainty  in 
this  measurement.  On  this  basis,  it  would  seem  that  las- 
ing or  Raman  gain  at  the  idler  could  make  a significant 
contribution.  However,  the  temporal  dependence  of  the 
signal  exhibits  little  of  the  spikey  or  noisy  appearance 
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associated  with  super  radiant  emission;  the  modulation 
of  the  signal  instead  follows  that  of  the  laser  quite  well 
as  expected  for  parametric  gain.  Additionally,  if  the 
signals  are  attributed  to  lasing  at  the  idler  frequency 
it  is  estimated  that  the  strength  of  the  idler  at  which 
saturation  due  to  depletion  of  9D^  atoms  should  occur  is 
about  500  watts  for  the  3 ym  line  instead  of  approximately 
.1  watt. 

A third  possible  lasing  process,  two-photon  stimula- 
ted emission  from  state  9D^  to  6S,  would  produce  signal 
and  idler  frequency  photons  because  of  the  resonant  en- 
hancement provided  by  the  P states.  This  process  is 
ruled  out  because  no  signal  wave  is  observed  in  the  re- 
verse direction  and  because  population  inversion  between 
9D^  and  6S  is  unlikely. 

We  conclude  that  lasing  and  Raman  gain  are  not  domi- 
nant in  the  production  of  the  radiation  we  have  studied. 
That  is,  the  contribution  of  population  inversions  to  Ak 
is  not  so  large  that  the  process  resembles  lasing  although 
population  shifts  are  large  enough  to  create  power  depen- 
dent changes  in  Ak. 

Our  analysis  has  been  based  on  a steady  state  model 
and  the  further  idealizations  of  a parallel  undepleted 
pump  beam,  no  significant  Stark  shifts,  no  diffraction 
losses  of  the  idler,  or  absorption  of  the  idler  by  CS2 
molecules,  or  ionization  losses.  In  principle  the  system 


could  be  described  completely  by  the  solution  of  the 
equations  of  motion  of  the  density  matrix  including 
these  effects  and  all  the  relevant  atomic  states  of 
cesium.  This  would  take  into  account  population  shifts 
generated  by  the  four  parametric  processes  as  well  as  lin 
ear  gain  or  loss,  Raman  processes,  ionization  losses, 
beam  geometries  and  so  on.  At  this  stage  we  are  content 
with  a qualitative  explanation  of  the  main  features. 

If  a more  detailed  investigation  is  desired,  it  should 
be  undertaken  in  a more  suitable  experimental  system 
in  which  competing  processes  and  complications  intro- 
duced by  loss  of  atoms  and  frequency  shifts  are  minimized 

5.4  Summary  & Conclusions 

We  have  made  the  first  detailed  measurements  of  the 
power  and  density  dependence  of  the  intense  resonant 
radiation  generated  in  an  atomic  vapor  pumped  by  a power- 
ful laser.  In  particular  we  have  studied  them  in  cesium 
vapor  pumped  by  a ruby  laser  timed  near  the  6S  - 9D| 
two-photon  resonance.  These  measurements  are  shown  to  be 
in  reasonable  agreement  with  a simple  steady  state  calcu- 
lation of  the  parametric  gain  for  small  signal  and  idler 
intensities.  Considering  the  extremely  complex  nature  of 
the  actual  interaction  which  involves  four  signal  and 
idler  fields  all  simultaneously  present  and  interacting 
with  one  another  through  the  various  state  populations. 
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and  the  possibility  of  lasing  between  various  levels 
and  possible  higher  order  processes,  these  calculations 
are  not  expected  to  yield  a precise  fit  to  the  data, 
but  we  believe  we  have  demonstrated  the  existence  of 
parametric  gain,  and  we  have  made  a detailed  investigation 
of  saturation  and  qualitatively  explained  the  saturation 
mechanism. 


APPENDIX  A 


STEADY  STATE  DENSITY  MATRIX  CALCULATION 

Expressions  describing  third  harmonic  generation 
and  two  photon  absorption  which  are  applicable  to  our 
experimental  situation  are  derived  in  this  appendix. 

The  density  matrix  formulation  is  used  and  a steady 
state  solution  is  assumed.  Doppler  averaging  is  in- 
cluded. The  model  quantum  system  upon  which  the  cal- 
culation is  based  consists  of  an  electromagnetic  field 
of  angular  frequency  o>  incident  upon  a multilevel 
quantum  system  consisting  of  levels  6S^,  9D and  the 

series  of  states  nP  (see  fig.  A.l) . Only  states  6S, 

*5 

and  9D3  are  allowed  to  have 
significant  populations. 
First  the  populations  of 
the  states  6S^  and  9D^  will 
be  calculated,  then  the  ex- 
pectation value  of  the  di- 
pole matrix  element  y(3w)  is 
derived.  The  former  is  re- 
lated to  the  amount  of  TPA; 
the  latter  to  the  THG  signal, 


nP 


' 

S ' / * ✓ 


9D 


A. 


6S)j 


fig.  A.l 


The  density  matrix  p is  defined  in  the  conventional  man- 
ner so  p is  the  population  in  state  n and  the  expecta- 
tion value  of  the  operator  A is  given  by  <A>  = Trace {pA} 


L29 
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TP  A 


The  equations  of  motion  used  to  find  p^2  (the  popu- 


lation in  state  9D3)  are 

'S 


ihp 


ih 


12 


hu)  P + fH,p"l  - m-  P 
12  12  L J12  T2  12 


ihp 


ih 


21 


”hwl2P21  + HH ' pH  2 1 “ T2  p2 1 


ihp 


11 


[H'P]1L  + T1  P22 


(A.l) 
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P11  + P22  1 


where  T.  is  the  lifetime  of  state  2,  T is  the  lifetimes 
1 ^ 


of  P12  and  P21'  and  H is  the  perturbation  Hamiltonian, 


H = 


-e 


12 
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I (E  e-iUt  + E*  eiwt) 2 


(A. 2) 


where  d = d and 
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If  a steady  state  solution  is  assumed  (Pi;L 
and  p has  the  form 


p22  °) 


-2itot 

■pll 

p12e 

p 2icot 

Lp2ie 

P22 

(A. 4) 


then  the  set  of  equations  A.l  may  be  solved  for  p.  The 
result  is 


-2  icot  . > 2 

1 9 d12  e (pll  “ p22 ^ E 


•(e/h) 


(A. 5) 


2w  - W]_2  + sr- 


P11  P22 


•J-l  4 1 12 

1 + sr-  (e/h)4  


d19|2  ( I E I 2 ) 2 


(2w  - c12)2  + (1/T2)2 


(A. 6) 


If  the  quantities  S and  D are  defined  as 


j 1^2  (e/h)  4 Id-,,12  ( | E | 2 ) 2 


D = (2co  - uj12)2+  (1/T2)2 


(A. 7) 


(A. 8) 


then  the  expression  for  (P-^  - P22)  and  p22  ma^  wr^t_ 
ten 


P11  ” P22 


(A. 9) 


1 + s/(dt2  ) 
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9 9 0 9 

* 2T2  D + S/T2 


(A. 10) 


When  2w  = u)  , the  expression  for  p becomes 
12  22 


22 


2 1 + S 


(A. 11) 


This  describes  the  population  in  state  2 for  exact 
resonance  tuning.  At  low  field  intensities  of  the  ap- 
plied field,  p22  is  proportional  to  the  square  of  the 
intensity  while  at  very  high  intensities  it  is  .5. 

When  the  incident  radiation  is  detuned  away  from  the 
resonance,  p^2  decreases  from  its  resonant  value.  Equa- 
tion A. 10  describes  the  shape  of  the  resonance.  It  has 
a Lorentzian  lineshape  with  the  FWHM  value  of  2/1  + S/T2 . 
Thus  the  linewidth  is  2/T 2 at  low  intensities  but  broad- 
ens in  proportion  to  the  intensity  when  S is  much  larger 
than  unity. 

Doppler  averaging  will  now  be  included.  The  Doppler 
averaged  form  of  eqn.  A. 10  is 


<P 


22 


> 


S _1 1_ 

/l  + S 2T2  Auj2 


z" (x, y) 


(A. 12 


where 


Au>2  E 2u)  { 2 kT/Mc^ } ^ 


(A. 13 


and  Z (x,y)  is  the  imaginary  part  of  the  Plasma  Dis- 
persion Function  with  the  arguments 


w12  - 


y = 


a + s 


T2  Aw2 


(A. 14) 


Finally,  note  that  the  fluorescence  from  state  2 is 
proportional  to  <p22>  so  t*ie  TPA  signal  is 


TPA  <*  <p  > = 

22 


2T2'Aw2 


(A. 15) 


The  polarization  at  the  frequency  3w  is  proportional 
to  the  3 co  component  of  the  expectation  value  of  the  di- 
pole transition  operator  y.  Since  the  expectation  value 
is  found  from 


<y>  = Trace  (py) 


(A. 16) 


we  must  find  the  third  harmonic  components  of  p.  This 
is  done  by  solving  the  equations  of  motion  of  p in  the 
presence  of  a second  perturbation  Hamiltonian  H2 . For 
simplicity  only  one  of  the  nP  states  is  considered.  It 
is  labeled  state  3.  Later  the  contributions  from  each 


of  the  nP  states  will  be  added  to  produce  the  total 
<y  ( 3gj)  > . 
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(A. 17) 


The  equations  of  motion  p are  displayed  below  where 
p , P22'  P 22  and  P21  are  assumed  to  have  the  values 
calculated  above. 
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These  are  solved  for  p31  and  p13  and  used  in 


<U(3uj)>  = Trace  (pp  ) = p p + p p (A. 19) 

13  13  31  31  13 


The  result  is 
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(A. 20) 
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The  terms  — ££_ 


(u  - 3w  + iT3_1 ) 


must  now  be  replaced  by  a 
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sum  of  similar  terms,  one  from  each  of  the  nP  states, 

plus  a continuum  contribution. 

Since  x ( 3w  ;to,  co  ,co)  E e^wt  is  the  second  term  of 

2 

eqn . A. 20  and  THG  « |<X>I  • 


THG  « I3  |<  5-  r >|2 

“ 1 + S/(DT0^)  - 2u  - iT 


12 


Performing  the  Doppler  averaging  yields9 


(A. 21) 


THG  - - | Z ' ^ (x ,y ) + -L-g  Z "2  (x,y)  } 
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where  D 


x,  y and  z are  defined  above. 


APPENDIX  B 


PHOTOIONIZATION  OF  9D3  AND  6S  STATES  IN  ATOMIC  CESIUM 

The  calculation  of  the  rates  of  photoionization  of 
the  9D3  state  of  cesium  by  14414  cm-'*'  radiation  and  of 
the  6S^  state  by  43243  cm  ^ radiation  (14414  cm'-*  is  the 
wave  number  of  the  ruby  laser  and  43243  is  its  third 
harmonic)  is  outlined  in  this  appendix. 

The  method  of  calculation  is  the  quantum  defect 
method  of  Peach.17  Equation  23  of  Peach  gives  the 
following  formula  for  the  photoionization  cross  section 
a of  an  atom: 


8aa^v  , 

, O . I 9 i 3 

c ~ ( £ ,V)  = — {l  + e \>  } 
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3z  C(v,£) 


(B.l) 


l C5  , [G  ( v £ ; e ' l ' ) Costt  { v+p  ' (e  ' ) + x(  v£. ; e ' £ ' ) }] 
' = £,  + ! * 


The  quantities  have  the  following  definitions: 

unprimed  state  is  bound  state;  primed  state  is 
ionized  state 

a = e^/ch  ; aQ  = .529><10-®cm 

v = frequency  of  radiation 

c'  = kinetic  energy  of  ejected  e in  Rydberg  units 

z = charge  of  core  = 1 for  neutral  cesium 

v = effective  principal  quantum  number 
2 2 

ie.  v = z /Inl 
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I U = threshold  ionization  energy  of  nJl  electron 
in  Rydbergs 

n = actual  principal  quantum  number 
y = n - y ; y is  called  the  quantum  defect 

y'(e')  is  the  extrapolated  quantum  defect 
in  the  continuum 

, = a coefficient  depending  on  l and  l' 


C (v, i) 
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V 

G and  X are  tabulated  by  Peach  for  l,  H < 3 

In  order  to  calculate  the  cross  section  for  a parti 
cular  case  one  calculates  v,  e',  and  C^,  and  estimates 
the  value  of  y'(e')  by  extrapolating  from  y(e)  (the 
quantum  defect  of  the  bound  states)  and  estimates  ^ (v , l) 
The  values  of  G(v£,;e'£.')  and  X(\>SLj  e'  H* ) are  interpolated 
from  the  tables  of  Peach.  The  results  are  expected  to 
be  correct  to  within  a few  percent  for  the  calculations 
of  interest  here. 

The  formula  above  did  not  include  fine  structure. 
This  can  be  included,  however,  by  using  the  following 
values  (which  are  correct  for  atoms  which  are  excited 
to  state  9D^  by  linearly  polarized  light  and  ionized  by 
light  of  the  same  polarization  or  for  atoms  ionized  from 
the  state  6S  by  polarized  light)  in  place  of  C0. 
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C(S^  -►  Continuum  P)  = 1 

C(D3  -►  Continuum  P)  = 17/25  (B.2) 

C(D3  -+  Continuum  F)  = 18/25 

The  values  of  the  ionization  rates  which  were  cal- 
culated in  this  manner  are: 


W (9D3  m.  = ±h  -*■  Continuum  P)  = 

h 3 

.457  I sec  ^ ; A = 693.5  nm 

W ( 9D3  m.  = Continuum  F)  = 

4 3 

7.16  I sec"^  ; A = 693.5  nm 


W (6S^  -+  Continuum  P)  = 

.33  I sec"^  ; A = 231.4  nm 
(I  in  watts/cm^) 


(B . 3) 
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APPENDIX  C 


A.C.  STARK  SHIFTS  OF  THE  6S.  AND  9D3  STATES 

h i 

OF  ATOMIC  CESIUM  BY  693.5  nm  RADIATION 


This  appendix  contains  the  outline  of  the  calcula- 
tion of  the  Stark  shifts  induced  by  the  ruby  laser  radia- 
tion field  (at  693.5  nm)  on  the  6S.  and  9D3  states  of 
cesium.  The  shift  of  a state  is  proportional  to  the 
linear  polarizability  of  that  state  so  there  is  a con- 
tribution to  the  shift  from  the  bound  electron  states 
and  from  the  continuum  states.  First,  the  bound  state 
contribution  is  calculated  for  both  states  (9D3  and  6S.  ) 
and  later  the  continuum  contribution  is  calculated. 
Linear  polarization  of  the  radiation  field  is  assumed 
and  only  the  shift  of  m^  = ±h  states  is  of  interest 
(states  of  different  m shift  by  different  amounts.) 

Bound  State  Contributions 

The  size  of  an  A.C.  Stark  shift  is  given  by 


“nJtJm.  * “nljmj  1 ' (1-05  * ltrll> 


(C.l) 


2 — 1 

where  I is  in  watts/cm  , AE  is  in  cm  , and  a is  in 

( A° ) ^ . a „ . will  be  written  a below  since  it  is  not 
nx- jm^ 

the  polarizability  as  usually  defined  but  is  the  polari- 
zability of  a particular  nr  state  for  linearly  polarized 
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light.  The  definition  of  a is 


2Ne‘ 


a = 


l { | <n  ' l ' s ' j 'm.  | z | nfl-s  jm . > | 
i 3 3 


(c.2: 


(u/)/(w'2-  w2)  } 


where  a)'  is  the  transition  frequency  between  the  initial 
and  final  states  and  the  sum  is  over  final  states. 

This  definition  is  to  be  contrasted  with  the  usual 
definition  of  polarizability  which  is  just  the  average 
of  5 over  the  m.  states,  ie.. 


01  “ 2j  + 1 E 01 

m . 

3 

It  can  be  shown  that  5 = Ca  where  C is  listed  be- 


(C . 3) 


low . 


Transition 


- Ph 


D 


D,  ^ 

2 "2 


c 

2 

1/5 

6/5 


Trans ition 


Q ->•  P . 


Now  a is  related  to  f numbers  by 


Ne 

a - — l 


f . 
x 


2 2 
m i u>  -u 
i 


C 

1 

1 


(C  . 4 ) 
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where 


j _ 8fT^mc  S 

3he2 A 2 j + 1 

S = £ | <n  ' 5- ' s ' j ' m ' | r ] n£s  jm  . > | ^ 

m . ^ 

1 

I 

m 

j 


(C  . 5) 


(C . 6 ) 


The  values  of  f for  the  various  transitions  can  be 
calculated  in  the  Coulomb  approximation  by  the  methods 
of  Bates  and  Damgaard1 8 or  taken  from  Stone.19  Then  a 
can  be  calculated  and  a follows . The  discrete  state 
contribution  to  the  ct's  calculated  in  this  manner  are 
listed  below. 


a9D, (m  =±h) 

h j 

oc6S 

*5 

a6S 

*5 


-28.43  A°3  ; 

Coul . 

f#'s, 

A = 693.5 

nm 

-100.02  A°3; 

Coul, 

f#'s, 

X = 693.5 

nm 

-116.03  A°3 ; 

Stone 

f#'s. 

A = 693.5 

nm 

Continuum  States  Contribution 


trom  the  Kramers  - Kronig  relations  it  can  easily 
be  shown  that 


a(uJo} 


oo 

_L 

X 4 7T  3 1097' 


a ( co  ) dco 


J0  4tt  j 109737  co ' 2 — co  2 
0 o 


(C . 7) 


j 3 

where  the  units  are  £0]  = cm  , Qo]  = Rydbergs , fa]  = cm  • 
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